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Abstract 


This  report-  introduces  a  new  method  for  solving  the  problem  of 
optimally  spacing  points  in  a  three-dimensional  region  so  that  their 
distances  from  each  other  are  as  great  as  possible.  One  application  of 
the  problem  deals  with  color  selection  for  aircraft  displays  where  the 
colors  are  plotted  as  points  in  a  three-dimensional  color  space  and  the 
distance  between  two  points  is  directly  related  to  the  di sti ngui shabi  1  ity 
of  the  two  colors.  The  method  itself  is  a  heuristic  algorithm  very 
similar  to  one  designed  by  Carter  and  Carter  (2).  The  newer  algorithm 
apparently  yields  similar  solutions  with  fewer  runs,  but  because  it  is 
more  thorough,  it  is  slower.  The  program  was  tested  on  problems  as  large 
as  23  points  whose  feasible  region  had  seven  faces.  The  major  dis¬ 
advantage  of  this  new  method  is  that  its  solutions  are  not  guaranteed  to 
be  optimal.  As  a  result,  the  user  must  perform  several  replications  of 
various  randomly  selected  starting  locations  in  order  to  increase  the 
chances  of  achieving  an  optimal  solution.  - 


I.  Introduction 


For  centuries,  mathematicians  have  tried  to  solve  tricky  little 
puzzles  that  perplex  and  confuse  the  human  mind.  Rubick's  Cube  is  a 
popular  recent  example  of  this  phenomenon.  This  thesis  is  another 
example.  Simply  put,  the  problem  addressed  herein  is  to  space  points  in 
a  three-dimensional  region  (like  a  box)  so  that  they  are  as  far  apart 
from  each  other  as  possible.  On  the  face  of  it,  the  task  seems  rela¬ 
tively  easy,  but  nothing  could  be  further  from  the  truth.  A  simple 
example  would  be  to  maximize  the  distance  between  just  two  points  in  a 
cube.  While  the  answer  to  that  problem  might  be  intuitively  obvious 
(place  the  points  in  opposite  corners),  add  one  more  point  and  the  solu¬ 
tion  is  not  so  apparent  and  requires  a  great  deal  of  mathematical  rigor. 

Faced  with  a  problem  like  this  in  real  life,  most  of  us  would  be 
content  to  settle  on  something  that  is  less  than  optimal.  Rut,  to  a 
pilot  who  needs  to  distinguish  between  friendly  and  hostile  forces  on  his 
aircraft  display,  the  importance  of  an  optimal  solution  takes  on  a  dif¬ 
ferent  meaning.  That's  why  this  project  is  so  useful,  and  I  sincerely 
hope  that  it  will  help  that  pilot  in  even  a  small  way. 

One  final  thought  before  moving  on.  Although  this  is  a  Master's 
thesis  on  what  is  potentially  a  very  technical  subject,  it  is  my  desire 
to  make  the  report  understandable  to  even  the  most  uninformed  of  laymen, 
and  perhaps  even  enjoyable  at  times.  I  do  not  believe  that  writings  at 
this  level  were  meant  to  be  boring  and  esoteric.  So,  if  you  cannot 
understand  it,  then  I  have  failed  and  1  should  be  shot  at  dawn.  Rut  if 
you  don't  try  to  understand  it,  then  you're  already  wearing  the  blindfold 
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1 1 .  Background 


Why  in  the  world  would  one  want  to  space  points  in  a  three- 
dimensional  region,  anyway?  Well,  there  are  many  reasons.  Say  you 
wanted  to  locate  speakers  for  a  quadraphonic  stereo  system  in  your  living 
room.  You  would  want  them  as  far  apart  from  each  other  as  possible 
because  you  don't  want  the  various  sounds  to  become  confounded.  Or  per¬ 
haps  you  own  a  warehouse  with  a  limited  number  of  security  cameras.  You 
would  probably  want  the  cameras  optimally  spaced  so  that  no  two  were 
filming  the  same  areas  of  the  warehouse.  How  would  you  go  about  solving 
this  problem  mathematically?  You  could  try  to  maximize  the  sum  of  the 
distances  between  the  objects,  but  that  might  result  in  a  solution  with 
an  unacceptably  small  distance.  A  more  desirable  solution  would  probably 
result  from  maximizing  the  distance  between  the  two  closest  objects. 
This  "maximin"  formulation  is  accepted  as  the  objective  for  the  spacing 
problem. 

The  spacif  ]  problem  is  very  similar  to  a  problem  known  in  the  world 
of  operations  research  as  the  obnoxious  facility  location  problem.  An 
example  of  an  obnoxious  facility  would  be  a  nuclear  waste  disposal  site 
which  is  most  desirable  when  it  is  located  far  away  from  cities.  Church 
and  Garfinxel  (5)  offer  a  solution  to  the  obnoxious  facility  location 
problem  where  the  facilities  can  be  placed  in  a  discrete  number  of  loca¬ 
tions  on  a  network,  but  the  differences  between  that  problem  and  the 
spacing  problem  preclude  successfully  adapting  Church  and  Garfinkel's 
method  to  the  spacing  problem.  First  of  all,  the  idea  is  to  locate 
obnoxious  facilities  as  far  away  as  possible  from  other  fixed  points. 
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whereas  the  spacing  problem  aims  to  locate  the  objects  as  far  away  from 
each  other  as  possible.  Second,  Church  and  Garfinkel's  facilities  can  be 
placed  in  only  a  finite  number  of  locations  on  a  network,  while  the 
objects  in  the  spacing  problem  can  take  on  an  infinite  number  of  loca¬ 
tions.  Finally,  the  spacing  problem's  objective  is  to  maximize  the  mini¬ 
mum  distance  between  points  (maximin),  as  opposed  to  Church  and 

■ 

Garfinkel's  objective  of  maximizing  the  median  distance  (maxian). 

Another  problem  similar  to  the  spacing  problem  is  the  location  prob-  j 

lem.  It  seeks  to  locate  objects  (like  warehouses)  as  close  to  other  ; 

fixed  facilities  (such  as  retail  stores)  as  possible.  Once  again,  much 
work  has  been  done  in  the  area  by  many  people,  most  notably  Charalambous 
(1  and  3),  Cooper  (6),  Juel  (8  and  11),  and  Love  (10,  11,  and  18),  but 

* 

any  applicability  to  the  spacing  problem  is  negligible.  As  before,  there 
is  the  problem  of  locating  facilities  in  relation  to  fixed  objects 
instead  of  each  other.  The  location  problem  is  also  concerned  with  mini¬ 
mizing  the  sum  of  the  distances  (minisum)  or  minimizing  the  maximum 
distance  (minimax),  not  with  maximizing  the  minimum  distance.  ' 

There  are  many  practical  applications  of  the  spacing  problem.  In  ; 

j 

addition  to  those  mentioned  before,  the  spacing  problem  could  be  applied  ■) 

toward  locating  MX  missile  silos,  spacing  mines  in  a  mine  field,  or  per-  j 

haps  placing  communications  satellites  in  space  for  maximum  coverage  of 
the  earth.  A  particularly  fascinating  application  of  the  spacing  problem 
is  in  the  area  of  color  spacing.  One  of  the  tasks  of  the  Air  Force  Aero-  j 

space  Medical  Research  Laboratory  (AFAMRL)  has  been  to  choose  colors  for 
aircraft  displays,  air  traffic  control  displays,  aeronautical  maps,  etc., 
so  that  all  the  colors  are  as  distinguishable  from  each  other  as  pos¬ 
sible.  When  one  is  dealing  with  a  color  Cathode-Ray  Tube  (CRT)  display. 
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the  CRT  s  red,  green,  and  blue  guns  are  used  additively  to  produce  vari¬ 
ous  colors.  Therefore,  any  color  can  be  uniquely  defined  by  three  param¬ 
eters.  They  are  the  luminances  (brightnesses)  of  the  CRT's  red,  green, 
and  blue  guns.  So,  each  color  produced  by  the  CRT  can  be  plotted  as  a 
point  in  three  dimensions  where  the  axes  represent  those  three  values.  A 
depiction  of  the  resulting  color  space  is  shown  in  Figure  1.  The  boun¬ 
daries  of  the  region  correspond  to  the  technological  limits  of  the  CRT. 
For  instance,  in  Figure  1,  the  maximum  luminances  are  50,  160,  and 
20  candellas  per  meter  squarred  for  the  red,  green,  and  blue  guns, 
respectively,  with  zero  being  the  minimum.  In  addition,  there  is  a  small 
region  near  the  origin  that  is  not  a  feasible  color  choice  for  this 
example  because  dark  colors  are  not  desired. 


Red 

* 


Figure  1.  Example  Feasible  Region  for  the  Color  Problem 


Unfortunately,  this  red,  green,  blue  color  space  is  not  per¬ 
ceptually  uniform."  That  is  to  say  that  the  perceived  difference  between 
the  red  and  the  purple  shown  in  Figure  1  is  not  the  same  as  the  differ¬ 
ence  between  yellow  and  white  even  though  their  distances  are  the  same. 
The  Commission  Internationale  de  I'Esirage  (CIE)  recommends  the  CIE 
L*u*v*  system  as  a  more  perceptually  uniform  color  space  (13).  The  L* 
axis  is  a  function  of  the  luminance  of  a  color,  while  the  u*v*  plane 
identifies  the  color's  position  in  a  transformation  of  the  standard 
u'v'  CIE  color  diagram  (see  Figure  2).  The  distinguishability  of  two 
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Figure  2.  Color  Chart  for  1976  CIE-UCS  Chromaticity  Diagram 


colors  in  the  CIE  l*u*v*  space  is  approximately  proportional  to  the 
Euclidean  distance  between  their  points.  Euclidean  distance  is  nothing 
more  than  our  common  sense  notion  of  distance  whose  formula  is  given  by: 


d(i.j)  =>](Li  -  L-)2  +  (u*  -  u])2  +  (vf  -  vj7 


(1) 


where 

d(i,j)  =  the  distance  from  point  i  to  point  j 

( L * ,  uf,  v* J  =  the  coordinates  of  point  i 

Carter  and  Carter  have  developed  a  computer  algorithm  to  solve  the 
color  spacing  problem  and  cite  several  military  applications  (2:2936). 
They  describe  how  the  method  could  be  used  to  choose  colors  for  strategic 
aircraft  displays  where  the  different  colors  correspond  to  friendly,  hos¬ 
tile,  and  neutral  forces  or  various  enemy  target  types.  As  another 
example,  an  air  traffic  controller's  display  can  show  airplanes  at  dif¬ 
ferent  altitudes  all  represented  by  various  colors.  Some  engineering 
applications  include  showing  "the  distribution  of  some  property 
throughout  a  system,  such  as  stress  in  a  structure  or  percent  of  full 
capacity  in  various  parts  of  an  electric  power  grid"  (2:2936).  As  you 
can  see,  there  are  many  interesting  applications  of  the  color  spacing 
problem  both  in  the  military  and  in  the  private  sector.  AFAMRL,  however, 
has  not  been  able  to  successfully  implement  Carter  and  Carter's  program 
because  of  various  bugs  and  program  deficiencies  and  has  given  only  edu¬ 
cated  answers  when  decisions  regarding  color  selection  were  required  in 
the  past.  Concerning  the  algorithm.  Carter  and  Carter  themselves  believe 
that  "presumably  a  more  efficient  one  could  be  devised  by  specialists  in 


7 


operations  research"  (2:2937).  For  those  reasons,  this  effort  was  under¬ 
taken.  The  genera’  objectives  herein  are  to  provide  AFAMRL  with  a 
working  program  that  produces  good  answers  to  the  color  spacing  problem, 
and  is  also  an  improvement  over  Carter  and  Carter's  method. 
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III.  Li terature  Review  of  Carter  and  Carter's  Method 


The  only  known  solution  to  the  color-spacing  problem  was  introduced 
by  Carter  and  Carter  in  1982  (2).  It  consists  of  first  randomly  placing 
n  points  in  the  region  (shown  in  Figure  1),  then  identifying  the  minimum 
CIE  L*u*v*  distance  between  all  n(n-l)/2  pairs  of  points  in  that 
region.  Let's  call  that  value  D.  Once  that  distance  has  been  identi¬ 
fied,  the  two  closest  points  { i . e - ,  the  endpoints  associated  with  the 
minimum  distance)  are  investigated  to  see  what  effect  is  created  by 
moving  each  endpoint  to  its  26  adjacent  locations.  This  step  can  con¬ 
ceptually  be  thought  of  as  "wiggling"  the  two  closest  points  to  see  if 
they  can  be  moved  farther  apart.  Figure  3  shows  the  various  alternatives 
for  one  endpoint.  The  alternatives  fall  on  the  boundaries  of  a  cube 
whose  sides  are  twice  the  step-size  in  length  and  whose  center  (the  27th 
point)  is  the  endpoint  itself.  There  are  52  different  moves  that  can  be 


Figure  3.  Alternative  Locations  for  an  Endpoint 
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made  (26  for  each  endpoint).  If  a  move  decreases  the  distance  between 
those  two  points  or  pushes  one  of  them  outside  of  the  region,  that 
alternative  is  no  longer  considered.  The  remaining  alternatives  are  then 
ranked  according  to  how  much  they  increase  the  distance  between  the 
endpoi nts. 

If  the  highest  ranking  alternative  causes  an  increase  in  D,  that 
move  is  made.  Otherwise,  the  second  highest  ranking  alternative  is  con¬ 
sidered,  and  so  on,  until  one  of  the  alternatives  causes  an  increase  in 
D.  If  none  of  the  alternative  moves  increase  0,  then  the  step-size  is 
halved  and  the  process  is  repeated;  otherwise,  the  point  is  moved  and  the 
process  is  repeated  until  an  expanding  move  cannot  be  made,  at  which  time 
the  step-size  is  halved. 

This  halving  continues  until  the  step-size  is  less  than  one  lumi¬ 
nance  unit  in  the  red,  green,  and  blue  color  space.  Once  the  step-size 
is  less  than  one  unit,  it  is  increased  back  to  its  original  value  and  the 
point-moving  and  step-size  halving  process  is  repeated,  "to  check  that 
the  value  of  D  arrived  at  is  not  merely  a  local  optimum"  (2:2938).  If 
the  points  remain  in  the  same  place  throughout  the  repetition,  then  "the 
...  configuration  of  points  is  assumed  to  represent  a  global  optimum 
value  of  0"  (2:2938). 

In  their  results,  the  authors  state  that  when  the  number  of  colors 
(n)  is  three,  the  solution  is  identical  for  each  random  placement  of 
points,  and  that,  in  general,  the  number  of  identical  solutions  decreases 
as  n  increases.  They  also  give  the  results  of  example  problems  where  3, 
4,  6,  10,  and  25  colors  were  spaced. 

Carter  and  Carter's  method  is  a  very  good  first  attempt.  Nonethe¬ 
less,  it  h  some  problems  as  well.  For  instance,  there  is  no  supporting 


proof  that  the  solutions  are  globally  optimal,  although  Carter  and  Carter 
assume  that  they  are  (2:2938).  In  the  article,  the  authors  summarize 
results  of  their  procedure  where  many  replications  are  made  for  placement 
of  3,  4,  6,  10  and  25  points.  For  the  case  of  six  points,  50  replica¬ 
tions  were  made  and  the  variance  among  the  values  of  D  arrived  at  was 
62.45,  where  the  maximum  value  of  D  for  all  replications  was  124.08 
(2:2938).  This  relatively  high  variance  points  to  the  fact  that  a  wide 
range  of  values  for  D  can  be  expected  for  any  one  replication,  and  that 
most  replications  are  not  close  to  even  the  maximum  known  value  of  0  for 
that  problem,  let  alone  the  global  optimum  value  which  for  color  spacing 
problems  is  largely  still  unknown.  Whether  or  not  the  solutions  repre¬ 
sent  even  local  optima  is  uncertain  as  well. 

Another  problem  is  that  there  is  no  validation  of  results.  For 
instance.  Carter  and  Carter  could  have  tested  the  algorithm  on  the  simple 
case  of  spacing  eight  points  in  a  cube  to  see  if  the  optimal  solution  (a 
point  in  each  corner)  is  obtained.  There  is  not  even  a  discussion  as  to 
the  physical  desirability  of  the  colors  that  were  generated  by  the 
algorithm. 

Also,  the  computer  code  is  poorly  documented  and  difficult  to  under¬ 
stand.  Variables  are  not  identified,  variable  names  are  often  conflict¬ 
ing  and  confusing,  and  statement  labels  are  poorly  numbered  and  sometimes 
not  found.  As  for  the  documentation,  the  number  of  internal  comments  is 
insufficient  to  make  the  program  easily  understandable. 

Carter  and  Carter's  method  has  its  good  points  as  well.  For 
example,  it  is  easy  to  understand.  Best  of  all,  the  technique  apparently 
yields  solutions  that  are  at  least  close  to  optimal,  which  is  all  they 
really  set  out  to  do  anyway.  But,  they  could  possibly  do  better.  For 


instance,  the  two  closest  points  may  be  optimally  spaced,  but  what  about 
the  rest  of  them?  Wouldn't  it  make  sense  to  go  through  the  algorithm 
once,  then  maximize  the  second  closest  distance,  the  third  closest,  and 
so  on,  until  all  the  points  are  optimally  spaced? 

Also,  Carter  and  Carter  maneuver  their  points  in  one  color  space  and 
then  transform  them  to  another  color  system  to  calculate  the  distances. 
I  believe  a  more  appropriate  approach  would  be  to  do  all  the  maneuvering 
and  calculations  in  the  same  coordinate  system.  This  would  not  only 
increase  the  efficiency  of  the  program,  it  would  increase  its  flexibility 
as  well.  It  would  enable  the  algorithm  to  solve  any  type  of  spacing 
problem,  not  just  the  color  spacing  problem.  An  algorithm  like  Carter 
and  Carter's,  but  with  better  validation  and  documentation,  including  the 
concept  of  successively  maximizing  the  minimum  distance,  second  minimum 
distance,  third  minimum,  and  so  on,  as  well  as  continuous  operation  in 
the  same  coordinate  system,  should  make  for  an  improved,  more  reliable, 
more  understandabl e,  more  efficient,  and  more  flexible  solution 
technique. 
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IV.  Mathematical  Development 


If  you're  afraid  of  mathematics,  like  so  many  of  us  are,  this  sec¬ 
tion  might  not  appeal  to  you  at  first.  But  don't  be  frightened,  because 
the  mathematics  here  are  really  nothing  more  than  algebra,  although  the 
terminology  and  notation  might  be  somewhat  difficult  to  keep  track  of. 
The  first  part  will  show  the  mathematical  formulation  of  the  color- 
spacing  problem  in  general  terms.  Next,  there  will  be  an  examination  of 
some  of  the  solutions  to  simple  spacing  problems.  The  final  section  will 
be  devoted  to  discussing  how  the  problem  could  be  solved  by  nonlinear 
programming. 

Problem  Formulation 

For  the  spacing  problem,  we  want  to  maximize  the  minimum  distance 
among  n  points  in  a  convex  polyhedron.  It's  a  mouthful,  isn't  it?  Let's 
take  it  apart  and  define  some  terms.  A  polyhedron  can  be  simply  defined 
as  a  region  bounded  by  "many  faces."  A  triangle  is  an  example  of  a  poly¬ 
hedron  that  has  three  faces;  a  square  has  four  faces.  In  three  dimen¬ 
sions,  a  cube  and  a  box  are  each  polyhedrons  with  six  faces  while  a 
pyramid  has  five.  These  are  all  convex  regions  as  well,  because  there 
are  no  "juts"  sticking  into  the  region.  A  star  (for  instance,  one  of 
those  on  the  U.S.  flag)  is  not  a  convex  region  for  that  very  reason. 
More  formally,  a  region  is  convex  if  every  point  in  that  region  can  be 
connected  to  every  other  point  in  the  region  with  a  straight  line  that 
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stays  completely  within  the  region.  All  it  means  is  that  no  juts  or 
dents  are  allowed. 


In  two  dimensions,  the  region  bounded  by  a  circle  is  convex,  as  is 
that  of  a  sphere  in  three  dimensions,  but  neither  one  of  them  is  a  poly¬ 
hedron  because  they  have  no  faces.  Mathematically  speaking,  a  polyhedron 
must  have  linear  boundaries.  That  means  that  a  polyhedron  must  be  made 
up  of  a  series  of  straight  lines  (or  planes  if  it's  three-dimensional). 
No  curved  boundaries  are  allowed.  Since  the  color  spacing  problem  is  a 
three-dimensional  one,  we  will  largely  concern  ourselves  with  three 
dimensions  from  now  on. 

The  boundaries  are  known  in  the  trade  as  constraints.  A  linear  con¬ 
straint  is  of  the  form: 


ax  +  by  +  cz  =  k 


(2) 


where 

(x,y,z)  =  the  coordinates  of  a  point 
k  =  a  constant  known  as  the  right  hand  side  value 
a,  b,  and  c  =  any  real  numbers 

The  above  equation  describes  a  plane.  Changing  the  equation  to  ax  + 
by  +  cz  k  ,  we  have  the  points  on  a  plane  in  addition  to  those  on  one 
side  of  it.  Rut  our  region  has  many  faces  to  it,  so  that  there  are 
actually  many  constraints  of  the  form: 

aix  +  biy  +  cl*  1  kl 
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a2x  +  b2y  +•  c2z  _<_  k2 

a3x  +  b^  +  C3Z  <_  k3 

•  *  •  • 

•  0  •  • 

•  •  •  • 

afx  +  bfy  +  cfz  _<  kf 

where 

f  =  the  number  of  faces  which  bound  the  region 

Together,  these  constraints  form  a  region  within  which  all  of  our  points 
must  stay.  If  a  point  is  outside  of  the  region,  then  at  least  one  of 
these  constraints  will  be  violated. 

Let's  turn  our  attention  again  to  the  objective.  We  want  to  maxi¬ 
mize  the  minimum  distance  among  n  points.  The  distance  from  one  point  to 
another  in  the  region  can  be  measured  by  the  Euclidean  distance  formula: 


where 

d ( i , j )  =  the  Euclidean  distance  from  point  i  to  point  j 
(xj,  y^,  z.j)  =  the  coordinates  of  point  i 

In  order  to  solve  the  problem,  we  need  to  find  the  minimum  distance 
between  all  pairs  of  points.  If  there  are  only  three  points,  it's  easy; 
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there  are  only  three  pairs  of  points  and  we  want  the  min  (d(l,2), 
d(l,3),  d ( 2 , 3 )  }  where  min  {  }  denotes  the  smallest  value  of  all  the 
numbers  in  the  set  in  brackets.  But  for  n  points,  that  set  may  be  quite 
large.  The  minimum  distance  for  n  points  would  be: 

min  {d(l  ,2) ,  d(l,3),  d(l,4),  ....  d(l,n),  d(2,3),  d(2,4),  ..., 

d(2,n),  d(3,4),  ...,  d(3,n),  ...,  d(n-2,n-l),  d(n-2,n),  d( n-1 ,n )  } 

As  a  matter  of  fact,  there  are  n{n-l)/2  values  in  the  set.  Letting 

D  =  min  {d(i ,j) },  the  objective  then  is  to: 

maximize  0  =  maximize  min  {d ( i , j )  }  (4) 


subject  to 

4"(x1_-  x2)2  +  (yt  -  y2)2  +  (zj  -  i2f  “  <1(1 ,2) 
>  (xj  -  x3)2  +  (yx  -  y3)2  +  (zj  -  z3)2  =  d(l ,3) 


^Vl  -  +  ^n-1  •  +  (zn-l  -  Zn^2  =  d(n'1’  n) 


a1x1  +  b1y1  +  c1z1  _<  kL 


alx2  +  ^2  +  clz2  —  kl 
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n  =  the  number  of  points 
f  =  the  number  of  faces 

the  first  set  of  constraints  =  the  distance  equations 
the  second  set  of  constraints  =  the  boundary  constraints  for 
all  n  points 

We  can  eliminate  the  square  root  signs  in  the  distance  constraints 
and  it  will  not  affect  the  formulation,  so  let's  do  that  just  to  uncom¬ 
plicate  the  equations.  Also,  we  can  subtract  0  from  all  the  distance 
constraints  to  create  equations  of  the  form: 

-  Xj  )2  +  [yi  ~  Yj  f  +  (zi  -  Zj  f  -  D  =  d(i  ,j)  -  D  (5) 

But  we  know  that  d(i,j)  -  0  is  going  to  be  a  number  greater  than  or 
equal  to  zero  because  D  represents  the  smallest  possible  value  of  any 
d(i,j).  Therefore,  we  can  change  the  equation  above  to: 

(xi  -  xj  f  +  (yi  -  yj  f  +  (zi  -  Zj  f  -  D  >  0  (6) 

The  problem  now  becomes: 

maximize  0  (7) 

subject  to 

I*,  -  Xj)  2+  (y i  -  yj)2  +  (zi  -  z.f  -  D  >_  0 

alxg  +  bl^g  +  clzg  —  kl 
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a2xg  *  Vg  *  c2zg  i  k2 


afxg  *  bf*g  +  cfzg  <  "f 

where 

i  =  1  to  n  -  1 
j  =  i  +  1  to  n 

g  =  a  subscript  for  the  1  through  n  points  that  must  stay  within 
each  boundary  constraint 
n  =  the  number  of  points 
f  =  the  number  of  faces 

At  this  point  you  might  think  we're  done,  but  there  are  a  few  more 
simplifications  to  go.  For  practically  all  commercial  methods  of  con¬ 
strained  optimization,  the  constraints  are  required  to  have  equal  signs 
instead  of  inequalities.  We  accomplish  this  by  using  variables  that  are 
commonly  referred  to  as  surplus  and  slack  variables.  For  example,  we 
know  from  before  that  the  value  of  d ( i , j )  -  D  is  non-negative,  but  we 
don't  know  how  large  it  is.  If  we  call  this  value  s^j  and  subtract  it 
from  the  distance  constraint,  we  have: 

(x.  -  x.  f  +  (y.  -  y.  )^  +  (z-  -z-]^"D-s..  =  Q  (8) 

^  i  j  '  1  '•i  j  ;  ij 
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The  constraint  is  now  an  equality  constraint  as  required,  and  s^j  is 
called  the  surplus  variable  because  it  represents  how  much  greater  than 
zero  equation  (6)  is. 

Similarly,  we  can  add  what  is  called  a  positive  slack  variable  to 
the  boundary  constraints  to  make  equations  of  the  form: 

ahxg  +  Vg  +  chzg  +  Shg  =  kh  (9) 


where 

Shg  =  the  slack  variable 
The  problem  can  then  be  formulated  as: 


maximize  D  (10) 


subject  to 

h  -  *j)2  ”  "  (zi  -  zP2  -  ”  •  sij  ■ 0 

ahxg  +  Vg  *  chzg  *  Shg  =  kh 
V  *g*  zg*  °»  sij’  Shg»  kh  i  0 


where 

i  =  1  to  n  -  1 
j  =  i  +  1  to  n 

g  =  a  subscript  for  the  1  through  n  points 
h  =  a  subscript  for  the  1  through  f  boundary  constraints 
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n  =  the  number  of  points 
f  =  the  number  of  faces 


s^j  =  surplus  variables  for  the  distance  constraints 
Shg  =  slack  variables  for  the  boundary  constraints 

The  third  set  of  constraints  states  that  x,  y,  z,  D,  s,  S,  and  k 
must  all  be  greater  than  or  equal  to  zero,  commonly  referred  to  in  the 
business  as  the  non-negativity  constraints.  This  is  another  requirement 
of  the  commercial  optimization  techniques.  Let's  talk  about  the  non¬ 
negativity  constraints.  We  already  know  that  D  is  positive  (you  can't 
have  a  negative  distance),  and  s^j  and  Shg  are  also  positive  as  we 
discussed  before.  If  the  x's,  y's,  and  z's  could  possibly  come  out  nega¬ 
tive,  we  could  make  transformations  of  the  form  x-  =  x^  -  t  where  t  is 
the  minimum  possible  value  of  x,  to  ensure  that  x."  is  positive.  If  k 
is  negative,  then  you  can  multiply  both  sides  of  the  equation  by  -1  to 
alleviate  the  problem,  but  the  sign  changes  from  less  than  to  greater 
than,  and  Shg  becomes  a  surplus  variable  instead  of  a  slack  variable. 
These  steps  ensure  that  S  would  always  be  positive  as  well  as  k,  and  the 
non-negativity  constraints  are  now  completely  satisfied. 

Notice  that  we  now  have  n(n-l)/2  distance  constraints  that  are 
nonlinear,  and  n  •  f  linear  boundary  constraints.  The  next  step  is  to 
incorporate  the  nonlinear  constraints  into  the  objective  function  with 
what  is  called  a  Lagrangian  function.  If  we  call  the  left  side  of  the 
first  distance  constraint  Q^,  the  second  Qjj,  and  so  on,  the  Lagrangian 
would  look  like  this: 

L  =  D  +  X12°I2  +  X13°13  +  +  Xln°ln+  ***  +  xn-lnVln 
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or  similarly: 


n-1  n 

L  =  D  +  E  E  xij  Kx-j 

i =1  j=i+l  J 


xj)2  +  Oi  -  yj)2  +  ”  -  *ij] 


(12) 


where 


L  =  the  Lagrangian  function 
A.jj  =  the  Lagrangian  variables 

So,  the  formulation  of  the  problem  now  looks  like  this: 

maximi ze 


n-1  n 

L  ■  o  +  E  E 

i=l  j=i+l 


ij 


[(«, 


Xj)2  *  (y,  -  yj)2  +  (z,  -  Zj)2  -  H  -  5ij] 


(13) 


subject  to 

ahxg  +  Vg  +  chzg  +  Shg  =  kh 
x,  y,  z,  D,  s,  S,  k  >  0 


where 

g  =  a  subscript  for  the  1  through  n  points 
h  s  a  subscript  for  the  1  through  f  boundary  constraints 
i  *  1  to  n  -  1 


j  =  i  +  1  to  n 
n  -  the  number  of  points 
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f  =  the  number  of  faces 


s^j  =  surplus  variables  for  the  distance  constraints 
Shg  =  slack  variables  for  the  boundary  constraints 
(Xg,  yg,  Zg)  =  coordinates  of  point  g 
A .j j  =  Lagrangian  variables 

And  that's  it.  We  now  have  a  nonlinear  objective  function  and  n  •  f 
linear  constraints.  Additionally,  the  objective  function  is  mathemati¬ 
cally  classified  as  convex  because  it  is  the  summation  of  many  convex 
functions.  This  is  somewhat  unfortunate  because  if  it  were  concave,  the 
problem  would  fall  under  a  category  of  problems  which  could  be  solved 
using  convex  programming;  a  solution  technique  which  is  widely  used  in 
operations  research.  Refore  we  leave  this  section,  I  would  like  to 
mention  one  more  thing.  Originally,  it  was  thought  that  the  color  spac¬ 
ing  problem  had  completely  linear  boundary  constraints  (i.e.,  that  the 
region  was  a  polyhedron).  It  turns  out  that  this  is  not  the  case.  It  is 
indeed  a  convex  polyhedron  in  the  red,  green,  blue  coordinate  system,  but 
in  the  one  that  counts,  the  CIE  l*u*v*  system  where  distances  are  com¬ 
puted,  it  is  not.  However,  it  is  still  convex  and  the  general  formula¬ 
tion  above  still  holds.  Those  constraints  that  are  nonlinear  can  be 
moved  into  the  Lagrangian  function  in  much  the  same  manner  as  the  dis¬ 
tance  constraints  were,  while  the  linear  boundary  constraints  can  stay 
put.  That  concludes  the  problem  formulation,  and  we  are  now  ready  to 
look  at  the  solutions  to  some  easy  spacing  problems. 
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Solutions  to  Simple  Problems 

This  is  the  section  where  the  puzzle  is  solved.  Although  the 
solutions  here  are  for  very  simple  problems  and  many  of  them  can  be 
recognized  intuitively  without  any  trouble,  some  of  them  are  not  so 
apparent.  What  I'm  going  to  show  you  here  are  the  solutions  for  spacing 
2,  3,  4,  5,  and  6  points  in  a  square  to  get  a  feeling  for  the  types  of 
solutions  that  can  be  expected  from  the  spacing  problem.  Included  is  a 
discussion  of  the  mathematical  aspects  of  some  of  the  solutions  and  some 
conclusions  based  on  what  we've  seen. 

Taking  a  look  first  at  the  solution  for  two  points  (shown  in 
Figure  4),  we  can  see  that  it  is  exactly  what  one  would  expect.  The 
points  are  in  opposite  corners.  Realize  also  that  there  are  actually  two 
possible  optimal  solutions  that  will  yield  the  same  value  for  D.  One  is 
pictured  here  and  the  other  would,  of  course,  have  the  points  in  the 


other  two  corners. 

The  solution  for  spacing  three  points  is  somewhat  less  obvious  (see 
Figure  5).  It  turns  out  that  the  points  form  an  equilateral  triangle 


Figure  4.  Solution  for  Spacing  Figure  5.  Solution  for  Spacing 

Two  Points  in  a  Square  Three  Points  in  a  Square 


(i.e.,  a  triangle  whose  sides  are  of  equal  length)  with  one  point  in  a 
corner  and  the  other  two  points  on  opposite  faces  at  angles  of  15  degrees 
from  the  corner  point.  And,  not  coincidentally,  their  coordinates  are 
(0,0),  (1,  tan  15  degrees)  and  (tan  15  degrees,  1),  respectively,  for  a 
square  whose  sides  are  one  unit  long.  In  addition,  one  can  easily  tell 
that  this  problem  actually  has  four  solutions  that  are  multiple  optimal; 
each  one  corresponding  to  a  different  corner  point. 

Spacing  four  and  five  points  in  a  square  is  trivial.  For  four,  the 
points  belong  in  the  corners.  The  five-point  solution  is  identical  to 
the  four-point  solution  with  the  additional  point  going  in  the  center  of 
the  region  (see  Figures  6  and  7).  Incidentally,  both  of  these  problems 
only  have  the  one  optimal  solution. 


'  \ 


Figure  6.  Solution  for  Spacing  Figure  7.  Solution  for  Spacing 

Four  Points  in  a  Square  Five  Points  in  a  Square 


Now  we  get  to  the  really  exotic  one.  The  optimal  solution  for 
spacing  six  points  in  a  square  is  shown  in  Figure  8.  Is  that  what  you 
would  have  expected  the  solution  to  look  like?  If  you're  like  most 
people,  probably  not.  I  would  have  guessed  that  all  four  corner  points 
would  be  filled  and  that  the  two  other  points  would  be  somewhere  in  the 

interior.  Instead,  there  are  two  corner  points,  one  interior  point,  and 

three  points  on  the  faces.  The  presence  of  a  counterintuitive  solution 
like  this  one  indicates  that  accepting  human  judgment  for  good  solutions 
may  be  risky.  It  further  emphasizes  the  need  for  a  computer  generated 
solution  technique  like  the  one  introduced  here.  Notice  that  the  left 
two-thirds  of  the  region  is  identical  to  a  skinny  five-point  solution, 
and  that  the  points  in  the  right  two-thirds  form  a  diamond  where  all  the 

sides  are  of  equal  length.  In  all,  there  are  six  distances  that  are 

identical — the  four  of  the  diamond  and  the  two  connecting  the  interior 
point  to  the  corner  points — and  they  all  correspond  to  the  minimum 
distance.  Notice  also  that  there  are  four  multiple  optimal  solutions  for 
this  proDlem. 


Figure  8.  Solution  for  Spacing  Six  Points  in  a  Square 


That  brings  us  to  the  conclusions  that  can  be  drawn  from  these 
simple  examples.  First,  it  appears  that  the  points  generally  tend  to 
fall  in  the  corners  and  the  faces  before  they  go  to  the  interior  of  the 
region.  Second,  the  solutions  are  not  always  what  one  would  expect.  The 
six-point  problem  yields  a  good  example  of  that  phenomenon.  Third,  most 
of  the  problems  have  multiple  optimal  solutions,  so  that  the  decision 
maker  can  often  choose  between  several  alternatives.  Finally,  it  appears 
that  the  optimal  solutions  all  have  a  high  percentage  of  their  distances 
identical,  and  equal  to  the  maximized  minimum  distance.  For  instance, 
the  three-point  solution  has  3  out  of  the  3  equal  to  D,  the  four-point 
solution  has  4  out  of  6,  the  five-point  has  4  out  of  10,  and  the  six- 
point  has  6  out  of  15.  This  means  that  one  might  be  able  to  tell  how 
good  his  or  her  solution  is  not  only  by  the  value  of  D,  but  also  by  the 
number  of  distances  that  are  equal  to  that  value. 

Solution  by  Nonlinear  Programming 

Up  to  now,  you  may  have  been  thinking  that  the  problem  looks  inter¬ 
esting,  but  the  mathematical  formulation  appears  pretty  complex.  How 
does  one  actually  solve  it?  Well,  there  are  several  ways.  You  can  write 
a  heuristic  algorithm  like  Carter  and  Carter  did,  and  like  I  did.  Or  you 
can  solve  it  intuitively  if  it's  simple  enough.  Or  you  might  be  able  to 
use  a  computerized  nonlinear  optimization  technique.  This  section  will 
discuss  one  such  nonlinear  programming  method. 

It  was  originally  called  the  Method  of  Approximation  Programming 
(MAP)  when  it  was  first  introduced  back  in  the  late  1950s  by  Griffith  and 
Stewart  (7),  hut  now  it  is  perhaps  better  known  as  successive  linear  pro¬ 
gramming  (SLP).  It  works  like  this:  Each  nonlinear  constraint  is 


linearized  using  a  first-order  Taylor  s  series  expansion.  For  those  of 
you  who  don't  know,  the  Taylor's  series  expansion  of  an  equation  is  just 
a  mathematical  series  of  summed  terms  which  approximates  the  original 
equation  and  converges  to  the  original  with  enough  terms  in  the  series. 
Once  the  nonlinear  constraints  are  linearized,  the  user  selects  a 
starting  set  of  values  and  solves  the  now  completely  linear  problem  using 
linear  programming  (LP).  The  solution  to  that  linear  programming  problem 
is  then  used  as  the  new  starting  set,  and  another  LP  solution  is 
generated.  This  goes  on  until  the  LP  solutions  become  identical,  at 
which  time  the  LP  solution  is  also  the  solution  to  the  nonlinear  problem. 

There  are  certain  restrictions  on  the  type  of  problem  for  which  this 
technique  will  guarantee  an  optimal  solution,  yet  Griffith  and  Stewart 

state  that  "problems  have  been  solved  with  MAP  which  do  not  fully  satisfy 
all  of  these  requirements"  (7:379).  The  reason  this  technique  was  not 
chosen  to  solve  the  color  spacing  problem  is  because  of  the  size  of  the 

problem.  For  example,  in  order  to  solve  Carter  and  Carter's  problem  for 

2b  points,  there  would  be  175  boundary  constraints  (some  linear,  some 
nonlinear),  300  nonlinear  distance  constraints,  and  76  decision  vari¬ 
ables.  It  is  uncertain  whether  an  LP  program  could  solve  that  problem  in 
a  reasonable  amount  of  time,  let  alone  solve  it  again  for  the  second 
closest  distance,  the  third  closest,  and  so  on,  as  the  heuristic  algo¬ 
rithm  introduced  in  this  thesis  is  able  to  do.  Although  this  problem  may 
not  seem  like  an  impossible  task  for  some  commercial  LPs,  the  primary 

objective  of  this  effort  is  to  provide  AFAMRL  with  a  technique  that 
works.  For  that  reason,  I  chose  to  solve  the  problem  with  a  heuristic 
algorithm  instead  of  SLP. 


V.  Description  of  Heuristic 


If  you're  reading  through  this  section  and  it  looks  vaguely  familiar 
to  you,  it's  not  a  case  of  precognition  on  your  part,  so  don't  call  the 
psychi atri st.  It's  because  the  heuristic  described  here  is  largely  pat¬ 
terned  after  that  of  Carter  and  Carter  which  was  described  in  Sec¬ 
tion  III.  However,  there  are  some  major  and  very  important  differences, 
so  don't  think  that  you  can  skip  to  the  next  section  because  you  already 
know  all  about  it.  On  the  contrary,  you  should  read  through  this  entire 
section  if  for  no  other  reason  than  to  please  me.  I  will  be  extremely 
pleased  if  you  read  the  first  portion  which  basically  describes  how  the 
algorithm  works  and  some  of  the  logic  behind  it.  And,  I  will  be  ecstatic 
if  you  read  the  second  portion  which  identifies  the  differences  between 
this  algorithm  and  Carter  and  Carter's. 

The  A1 gori thm 

Recall  that  the  purpose  of  this  thesis  is  to  maximize  the  minimum 
distance  among  n  points  in  a  convex  region,  and  then  as  a  kicker  to  suc¬ 
cessively  maximize  the  second  minimum  distance,  the  third  minimum,  and  so 
on,  until  all  points  are  optimally  spaced.  The  first  step  in  achieving 
this  objective  is  to  randomly  place  the  points  uniformly  throughout  the 
region.  For  the  color  spacing  problem,  it  would  look  like  the  one  shown 
in  Figure  9,  Recall  from  before  that  this  region  is  not  a  polyhedron, 
although  Figure  9  depicts  it  as  one.  The  actual  color  spacing  region 
would  have  many  of  its  boundaries  bowed  out  slightly.  Because  they  are 
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Figure  9.  Approximate  Feasible  Region  for  the  Color  Spacing 
Problem  in  the  C IE  L*u*v*  Color  System 


bowed  out  and  not  in,  the  region  is  still  convex,  so  the  solution 
technique  described  here  still  applies. 

In  step  2,  the  program  finds  the  minimum  distance  among  all  n(n-l)/2 
pairs  of  points,  and  the  endpoints  that  correspond  to  that  distance.  The 
third  step  is  to  define  the  27  alternative  locations  for  each  of  the  two 
endpoints  found  in  step  2.  These  alternatives  lie  on  a  unit  cube  around 
the  endpoint  whose  sides  are  twice  the  step-size  in  length.  Figure  3  on 
page  9  depicts  these  alternative  endpoint  locations. 

In  the  fourth  step,  the  alternatives  are  checked  to  see  if  an 
improvement  can  be  made  by  moving  the  endpoints.  The  program  does  this 


by  first  checking  the  alternatives  to  see  if  they  are  in  the  region. 
Then,  taking  one  legal  alternative  at  a  time,  it  checks  the  distance 
between  endpoints  to  see  if  that  value  is  greater  than  before.  Next,  it 
checks  the  new  overall  minimum  distance  to  see  if  that  is  greater.  If 
all  of  those  conditions  are  successfully  met,  then  the  endpoints  are 
moved  and  the  routine  checks  the  next  alternative  for  an  even  greater 
possible  improvement.  It  does  this  for  every  possible  combination  of 
alternative  endpoint  locations;  a  total  of  729  combinations.  If  several 
alternatives  have  the  same  overall  minimum  distance,  which  is  likely,  the 
tie  goes  to  the  alternative  producing  the  greatest  distance  between 
endpoints. 

The  program  then  returns  to  step  2  and  finds  the  new  overall  minimum 
distance  and  the  endpoints  that  correspond  to  that  distance.  This  loop 
continues  until  no  improvement  can  be  made  with  the  given  step-size.  At 
that  time,  the  step-size  is  halved  and  the  program  begins  looping  again 
starting  at  step  2.  This  halving  process  continues  until  the  step-size 
reaches  some  minimum  value  chosen  by  the  user  and  we  enter  step  5. 

The  fifth  step  controls  the  fixing  of  points  already  maximized  and 
the  status  of  the  program.  When  the  step-size  reaches  its  minimum  toler¬ 
ance,  one  of  the  endpoints  is  fixed  in  its  position,  the  step-size  is 
returned  to  its  original  value,  and  the  program  returns  to  step  2.  The 
only  difference  in  the  logic  now  (other  than  a  point  being  fixed)  is  that 
the  minimum  distance  calculated  in  step  2  can't  have  two  fixed  points  as 


its  endpoints.  This  enables  the  program  to  concentrate  on  the  second 
shortest  distance.  The  fixing  of  points  continues  until  all  the  points 
become  fixed,  at  which  time  they  are  all  unfixed  and  the  routine  starts 
all  over  again  at  step  2  as  if  the  current  location  scheme  were  the 


starting  positions.  If  the  program  runs  all  the  way  through  again  and 
the  points  haven't  moved,  then  the  program  is  done.  Otherwise,  the  pro¬ 
gram  keeps  going  until  there  is  no  change  in  location  for  one  complete 
iteration.  Figure  10  contains  a  flowchart  of  the  program's  logic,  while 
Appendix  C  contains  an  example  of  how  the  program  optimally  spaced  four 
points  in  a  square. 

Differences  with  Carter  and  Carter's  Method 

As  I  said  at  the  top  of  this  section,  there  is  a  striking  resem¬ 
blance  between  Carter  and  Carter's  method  and  the  one  described  here. 
However,  there  are  four  major  differences  which  I  would  like  to  discuss 
at  this  time. 

The  first  of  these  differences  lies  in  the  checking  of  alternatives 
to  see  what  expanding  move  should  be  made.  Carter  and  Carter  move  only 
one  endpoint  at  a  time  and,  therefore,  choose  between  only  52  possible 
endpoint  location  schemes  (26  for  each  endpoint).  On  the  other  hand,  the 
Roley  algorithm  moves  both  endpoints  at  once  and,  as  a  result,  has  a 
total  of  729  possible  endpoint  location  schemes.  It's  just  a  more 
exhaustive  search  for  the  best  possible  move,  that's  all. 

Along  the  same  lines,  Carter  and  Carter's  method  ranks  the  alterna¬ 
tives  as  to  how  much  they  increase  the  distance  between  endpoints  and 
then  chooses  the  highest  ranking  alternative  that  also  increases  the 
overall  minimum  distance.  On  the  other  hand,  my  program  chooses  the 
alternative  that  most  increases  the  overall  minimum  distance  with  ties 
going  to  the  alternative  which  causes  the  greatest  increase  in  the  end- 
point  distance.  The  difference  is  that  Carter  and  Carter's  algorithm  may 


Figure  10.  Flowchart  of  Roley  Heuristic 


not  choose  the  best  alternative  for  increasing  the  overall  minimum  dis¬ 
tance,  whereas  mine  does. 

Another  difference  has  to  do  with  the  thoroughness  of  the  spacing 
method.  Carter  and  Carter  maximize  the  minimum  distance  but  leave  the 
second  minimum,  third  minimum,  and  all  the  others  where  they  are.  My 
algorithm  not  only  maximizes  the  minimum  distance,  it  successively  maxi¬ 
mizes  the  second  minimum  distance,  the  third  minimum,  and  so  on  until  all 
points  are  optimally  spaced.  The  result  is  a  complete  optimization  of 
all  distances,  not  just  the  minimum  distance.  The  price  is  a  more  com¬ 
plex  and  slower  computer  algorithm. 

The  fourth  and  final  difference  has  to  do  with  the  respective  coor¬ 
dinate  systems  that  are  used  as  the  primary  operating  spaces.  Carter  and 
Carter's  technique  spaces  points  in  the  red,  green,  and  blue  color  sys¬ 
tem,  and  converts  the  colors  to  C IE  L*u*v*  color  coordinates  in  order  to 
calculate  spatial  distances.  My  program  acts  in  a  reverse  manner.  It 
maneuvers  points  and  calculates  distances  in  the  CIE  L*u*v*  system,  but 
checks  color  locations  in  the  red,  green,  and  blue  system  to  make  sure 
thay  are  within  the  boundaries  of  the  region.  It  was  originally  hoped 
that  all  operations  and  calculations  could  be  done  in  the  CIE  t*u*v* 
system  for  increased  efficiency,  but  the  equations  for  the  boundary  con¬ 
straints  could  not  be  derived  in  the  L*u*v*  system,  preventing  its  use  as 
a  location  checking  coordinate  system. 


VI.  Results 


Now  that  you  know  how  the  thing  works,  let's  find  out  how  well  it 
works.  There  are  several  questions  that  need  to  be  answered  in  this 
section.  Does  the  algorithm  guarantee  an  optimal  solution?  Does  it  work 
for  very  large  problems?  How  well  does  it  perform  in  comparison  to 
Carter  and  Carter's  method?  What  are  some  of  the  factors  to  which  the 
program  is  sensitive?  What  are  its  biggest  advantages?  What  are  its 
biggest  disadvantages?  These  are  the  questions  that  will  indeed  be 
answered  in  this  section. 


Optimality 

The  algorithm  does  not  guarantee  an  optimal  solution.  However,  it 
was  tested  on  a  number  of  problems  whose  solutions  are  known  and  the 
results  are  encouraging.  These  problems  include  spacing  three,  four, 
five,  and  six  points  in  a  square  and  spacing  eight  points  in  a  cube. 
Each  problem  was  tested  with  three  different  sets  of  randomly  placed 
points  and  the  algorithm  produced  the  optimal  solution  in  14  of  the  15 
problems.  That's  a  success  rate  of  93  percent.  The  only  unsuccessful 
attempt  involved  one  of  the  tries  at  spacing  five  points  in  a  square.  In 
that  instance,  the  minimum  distance  arrived  at  was  .638  units  compared 
with  the  optimal  solution  of  .707,  so  that  it  was  within  90  percent  of 
optimal.  In  general,  the  more  replications  that  are  performed,  the 
better  chance  of  success. 


Size  Limitations 


This  may  be  the  largest  downfall  of  my  program.  I  had  originally 
intended  for  it  to  accommodate  up  to  50  points  and  10  faces.  However, 
the  largest  problem  it  has  been  able  to  successfully  solve  was  the  color 
spacing  problem  with  23  colors  and  7  faces.  Larger  problems  exceeded  the 
time  limit  of  1000  CPU  seconds  on  the  CDC  845  (Cyber)  computer.  The 
reason  for  this  size  deficiency  boils  down  to  the  thoroughness  of  the 
program.  Each  iteration  checks  n(n-l)/2  distances  for  up  to  729  alterna¬ 
tives.  It  takes  21  iterations  for  the  step-size  to  reach  its  minimum 
value.  The  program  must  perform  those  21  iterations  for  each  of  the  n 
points  that  are  fixed  in  succession.  It  then  goes  through  the  whole  pro¬ 
cess  at  least  one  more  time  or  until  the  point  locations  do  not  change. 
Experience  has  shown  that  it  usually  goes  through  twice.  With  an 
increase  in  n,  the  number  of  computations  (and  thus  the  CPU  time 
required)  goes  up  exponentially,  so  one  can  easily  tell  that  size  is  a 
serious  factor.  Presumably,  reducing  the  number  of  alternatives, 
increasing  the  minimum  step-size,  or  reducing  the  amount  of  point  fixing 
would  increase  the  ability  of  the  program  to  handle  larger  problems. 

Performance  i n  Compari son  to  Carter  and  Carter1 s  Method 

The  comparison  of  my  program  with  Carter  and  Carter's  was  done  with 
identical  color  parameters  of  Yg,  Ug,  Vg  ,  and  the  chromaticity  coor¬ 
dinates  of  the  guns.  Table  I  contains  comparative  results  for  spacing  3, 
4,  6,  10,  and  25  colors. 
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II YWIXi 


Number 

of 

Colors 

Best  Answer 
Method  (CIE  L*u*v*) 

Variance 

Number 
of  Repli 
cations 

Average  Computer 
Time  Per  Replication 

-  CDC  845 
(sec) 

UNIVAC  1182 
(sec) 

3 

Roley 

229.26 

1.17 

5 

23.50 

Carter 

239.33 

.11 

49 

1.31 

4 

Roley 

156.28 

136.62 

5 

42.95 

Carter 

155.87 

37.14 

50 

1.76 

6 

Rol  ey 

124.85 

79.61 

5 

100.49 

Carter 

124.08 

62.45 

50 

2.48 

10 

Roley 

71.07 

48.61 

20 

178.90 

Carter 

89.43 

27.73 

50 

4.22 

25 

Roley 

0.00 

0.00 

5 

>1000.00 

Carter 

51.60 

7.34 

46 

12.39 

As  an  explanation,  I  originally  intended  to  perform  only  five 
replications  for  each  problem  because  time  considerations  prevented  me 
from  performing  50  runs  as  Carter  and  Carter  did.  However,  I  ended  up 
doing  20  runs  of  the  ten  color  problem  to  try  to  decrease  the  great 
disparity  of  results.  Also,  I  was  unable  to  successfully  execute  the  25 
color  problem  because  the  CPU  limit  of  1000  seconds  was  exceeded  for  each 
attempt. 

It  is  difficult  to  compare  this  algorithm  with  Carter  and  Carter's 
because  there  are  so  many  factors  that  can  be  used  as  a  measure  of  effec¬ 
tiveness.  Most  people,  however,  would  consider  the  minimum  distances 
produced  by  each  method  as  their  primary  concern.  Notice  in  Table  I  that 


the  minimum  distances  are  very  similar  for  spacing  3,  4,  and  6  colors 
although  Carter  and  Carter's  results  are  based  on  approximately  10  times 
more  runs.  The  case  of  spacing  10  colors  is  perplexing.  One  would 
assume  that  the  results  would  again  be  equivalent  because  the  techniques 
are  so  similar.  Apparently,  that  is  not  true.  My  heart  tells  me  that 
something  is  wrong,  but  it  is  difficult  to  confirm  Carter  and  Carter's 
results  because  those  actual  color  locations  were  not  published  and  are 
not  available  at  this  time.  Perhaps  an  even  greater  number  of  replica¬ 
tions  would  yield  a  better  solution.  In  the  meantime,  I  can  only  admit 
that  my  algorithm  is  deficient  for  spacing  larger  numbers  of  colors. 

Another  criteria  one  might  use  to  judge  the  two  methods  would  be  to 
compare  the  distances  between  the  second  closest  points,  third  closest, 
and  so  on.  This  would  indeed  be  an  appropriate  measure  because  the  fun¬ 
damental  conceptual  difference  between  the  two  methods  is  that  mine 
concentrates  on  successively  maximizing  all  distances,  whereas  Carter  and 
Carter's  maximizes  only  the  minimum  distance.  Presumably,  these  efforts 
should  have  paid  off,  but  it  is  difficult  to  judge  because  of  a  limited 
sample  size.  Table  II  presents  a  comparison  of  distances  for  the  six- 
color  spacing  problem.  I  would  have  liked  to  compare  other  problems  in 
addition  to  the  six-color  one  but  it  was  the  only  one  that  was  published 
and  available,  so  the  sample  size  isn't  quite  what  one  would  need  to  make 
a  definitive  judgment.  Notice,  however,  that  the  first  three  distances 
are  slightly  better  for  my  algorithm,  the  next  six  are  fairly  even,  while 
the  last  six  are  clearly  in  Career's  favor.  This  could  indicate  that, 
while  maximizing  the  second  and  third  minimum  distances  produces  better 
results  for  those  values,  it  may  cost  in  the  long  run  in  the  form  of 
smaller  values  for  the  greater  distances.  Notice  also  that  the  first 
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TARLE  II 


COMPARISON  OF  ALL  DISTANCES  FOR  THE  SIX-COLOR  PROBLEM 


ith 

Minimum 

Di stance 

Roley 

Carter 

ith 

Minimum 

Di stance 

Roley 

Carter 

ith 

Minimum 

Distance 

Roley 

Carter 

1st 

125* 

124 

6th 

129 

130* 

11th 

182 

218* 

2nd 

125* 

124 

7th 

136* 

133 

12th 

201 

237* 

3rd 

125* 

124 

8th 

138 

140* 

13th 

220 

239* 

4th 

125 

125 

9th 

149* 

147 

14th 

239 

249* 

5th 

125 

125 

10th 

171 

210* 

15th 

241 

263* 

♦Denotes  the  superior  value. 

through  sixth  distances  are  basically  the  same  regardless  of  whether  or 
not  they  are  maximized.  They  seem  to  have  been  automatically  maximized 
by  maximizing  the  minimum  distance.  That  indicates  that  successively 
maximizing  higher  distances  may  not  be  worth  the  extra  computer  time 
required.  These  results  also  support  the  observation  that  was  made  in 
Section  IV  that  optimal  solutions  will  generally  have  several  distances 
equal  to  the  minimum  distance.  It  is  a  characteristic  that  could  be 
extremely  useful  in  determining  whether  or  not  a  particular  solution  is 
close  to  optimal . 

One  might  also  be  impressed  with  my  method  if  it  could  take  Carter 
and  Carter's  best  solution  for  a  particular  problem,  use  it  as  a  starting 
location  scheme,  and  improve  upon  it.  Using  Carter  and  Carter's  six- 
color  solution  (because  it  was  the  only  one  available)  as  an  initial  set 
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of  colors  resulted  in  a  mild  improvement  of  the  minimum  distance  from 
124.026  CIE  L*u*v*  units  to  124.145  units.  Research  by  Carter  and  Carter 
has  shown  that  distinguishabil ity  between  colors  "deteriorates  rapidly 
when  the  distance  between  colors  is  about  40  CIE  L*u*v*  units,"  (2:2937) 
so  this  slight  improvement  of  .12  unit  is  inconsequential.  Perhaps  more 
importantly,  my  algorithm  improved  on  11  out  of  the  14  ''ther  distances 
with  one  unchanged  and  two  decreasing  slightly.  You  may  have  noticed 
that  the  original  distance  of  124.026  does  not  correspond  to  Carter  and 
Carter's  calculated  value  of  124.08  as  shown  in  Table  I.  Taking  their 
same  set  of  color  locations  in  red,  green,  and  blue  coordinates,  I  cal¬ 
culated  a  minimum  distance  of  124.026  in  L*u*v*  coordinates  instead  of 
124.08  (the  exact  L*u*v*  coordinates  were  not  published).  The  difference 
is  probably  due  to  roundoff  error.  8ut,  how  does  one  explain  how  I  could 
improve  on  their  optimal  solution  with  basically  the  same  technique? 
Accuracy.  Because  Carter  and  Carter's  minimum  step-size  is  one  unit, 
they  are  not  able  to  attain  the  same  level  of  accuracy  as  my  heuristic 
does  with  a  minimum  step-size  of  .0001  units.  Consequently,  the  Carters 
are  not  able  to  "snuggle"  their  points  into  the  best  locations.  All  this 
indicates  is  that  the  minimum  step-size  is  a  control  over  the  level  of 
accuracy  one  wishes  to  achieve. 

As  was  briefly  discussed  before,  the  question  of  computer  time  may 
also  be  an  important  performance  characteristic  by  which  to  judge  the  two 
methods.  Theoretically,  Carter  and  Carter's  method  should  be  consider¬ 
ably  faster  because  it  maximizes  only  the  minimum  distance  and  isn't 
concerned  with  the  second  minimum  distance,  third  minimum,  and  so  on. 
Looking  again  at  Table  I,  we  can  see  that  this  is,  in  fact,  the  case.  It 
should  be  remembered,  however,  that  comparing  run  times  on  two  different 


40 


computers  is  somewhat  like  comparing  apples  and  oranges  and  should  be 
taken  with  a  grain  of  salt. 


As  far  as  those  at  AFAMRL  and  I  are  concerned,  the  chief  criterion 
for  comparison  is  which  program  works.  I  personally  know  of  no  place 
where  Carter  and  Carter's  algorithm  is  currently  running  and  several 
places  where  people  have  tried  to  implement  it  and  have  not  been  able  to, 
AFAMRL  included.  On  the  other  hand,  my  program  is  one  for  one  so  far  so 
it  must  be  given  the  edge.  Interested  parties  will  be  delighted  to  find 
out  that  my  program  is  heavily  documented  internally,  easy  to  understand, 
utilizing  classic  structured  programming  concepts  and  written  in 
FORTRAN  77.  These  characteristics  give  my  program  the  added  advantage  of 
being  easier  to  modify  and/or  upgrade  by  others  who  may  be  so  inclined. 
A  listing  of  the  code  is  contained  in  Appendix  8  for  all  to  judge  for 
themsel ves. 

Briefly  summarizing  this  section,  we  found  that  Carter  and  Carter's 
method  works  equally  as  well  for  finding  the  maximum  minimum  distance, 
but  requires  more  runs.  For  the  six-color  spacing  problem,  theirs 
yielded  slightly  worse  values  for  the  second  and  third  minimums,  and  much 
better  results  for  the  greater  distances.  Carter  and  Carter's  program 
takes  less  time,  but  has  a  slightly  lower  degree  of  accuracy.  Unlike 
Carter  and  Carter's,  my  program  works,  is  well  documented,  and  easy  to 
understand. 

Sensi t i vi ty  to  Parameters 

The  program  introduced  in  this  report  is  highly  sensitive  to  four 
parameters.  They  are  the  maximum  step-size,  the  step-size  reducing 
scheme,  the  minimum  step-size,  and  the  initial  random  positioning.  The 


maximum  step-size  is  important  because  some  solutions  are  not  obtainable 
for  certain  step-sizes  and  point  locations.  Take,  for  instance,  the  case 
of  spacing  eight  points  in  a  cube  whose  sides  are  one  unit  long.  Suppose 
seven  of  the  points  are  in  the  corners  and  the  eighth  is  in  r.ne  very 
middle  of  the  cube.  The  eighth  point  should  migrate  to  the  empty  corner, 
but  unless  the  step  size  is  greater  than  one-third  of  a  unit,  the  minimum 
distance  will  decrease  if  it  tries  to  move  in  any  direction.  Essen¬ 
tially,  the  point  is  "locked  in"  to  its  location.  To  avoid  this  problem, 
I  recommend  using  a  maximum  step-size  equal  to  approximately  half  the 
distance  between  the  two  farthest  points  in  the  region.  This  way,  a 

point  located  in  the  middle  can  move  to  any  other  spot  in  the  region. 

Granted,  a  point  doesn't  necessarily  have  to  be  in  the  middle  to  be 
locked  in.  Let's  say  the  eighth  point  in  our  example  above  is  not 

located  in  the  exact  middle.  Let's  say  instead  that  it  is  located 
towards  the  empty  corner  but  less  than  one-third  of  a  unit  from  the 
middle.  A  step-size  of  .5  now  puts  the  point  out  of  the  region,  so  it 
still  can't  move  closer  to  where  it's  supposed  to  be.  If  the  step-size 
is  halved,  the  point  still  may  not  be  able  to  go  towards  the  corner. 
But,  some  other  step-size  reduction  scheme  like  reducing  it  by  four- 
fifths  might  send  the  point  to  the  optimal  solution.  In  that  sense,  the 
reduction  scheme  is  another  factor  which  could  affect  the  solution.  The 
problem  with  that  alternative  is  that  it  takes  quite  a  few  more  itera¬ 
tions.  For  instance,  to  go  from  .5  to  .0001  by  halves,  it  takes  13 

iterations.  To  go  from  .5  to  .0001  by  four-fifths  takes  39  iterations. 
More  iterations  means  more  computation  time  and,  thus,  a  less  desirable 
program.  The  Carters  and  I  both  used  a  halving  scheme,  but  I  would 
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recommend  that  any  user  experiment  with  this  parameter  to  find  out  what 
is  best  for  his  or  her  application. 

The  minimum  step-size,  although  a  factor  to  the  solution,  is  not  as 
important  as  the  two  parameters  discussed  above.  The  minimum  step-size 
controls  the  accuracy  of  the  solution.  For  three  decimal  places  of 
accuracy,  .0001  should  be  the  minimum  step-size.  For  two  places,  use 
.001,  and  so  on.  Increasing  the  minimum  step-size  can  cause  a  surprising 
decrease  in  the  number  of  iterations  required.  For  instance,  it  takes 
only  six  iterations  to  go  from  .5  to  .01  by  halves  as  opposed  to  13  iter¬ 
ations  to  go  to  .0001. 

The  final  parameter  that  significantly  affects  the  solution  is  the 
seed.  This  is  the  starting  value  that  is  used  in  the  random  number 
generating  sequence  which  calculates  the  starting  locations  for  the 
points.  One  can  intuitively  realize  the  importance  of  the  seed  by  under¬ 
standing  that  some  initial  location  sets  are  better  than  others.  A  bad 
location  set  is  shown  in  Figure  11.  None  of  the  points  in  that  example 


Figure  11.  Example  of  Locally  Optimal  Solution 


43 


can  move  anywhere  without  getting  closer  to  another  point,  and  yet  the 
solution  is  not  optimal.  It  is  called  a  "locally"  optimal  solution  as 
opposed  to  the  "globally"  optimal  solution  which  has  a  point  in  each  cor¬ 
ner.  Spacing  problems  apparently  tend  to  have  many  local  optima.  Chang¬ 
ing  the  seed  and  doing  numerous  replications  will  increase  the  chances  of 
finding  a  global  optimum. 

Advantages 

The  heuristic  algorithm  described  here  has  numerous  advantages. 

Here  is  a  list  of  some  of  them: 

1.  Flexibility.  The  program  works  for  any  spacing  problem 
including  the  color  spacing  problem.  It  has  a  myriad  of  uses 
in  the  field  of  optics,  as  described  by  Carter  and  Carter 
(2:2939).  It  can  even  be  extended  to  four  dimensions  for 
possible  use  in  the  field  of  physics.  It  can  accommodate 
distance  formulas  other  than  the  Euclidean  distance  formula  or 
convex  regions  besides  polyhedrons.  Its  uses  are  limited  only 
by  the  limits  of  the  imagination. 

2.  Simplicity.  It  is  simple  to  understand  and  simple  to  use.  The 
structured  FORTRAN  programming  format  makes  it  simple  to  modify 
or  debug. 

3.  Reliability.  It  has  been  tested  favorably  against  simple 


problems  with  known  solutions  and  more  complex  problems  like 
Carter's  color  spacing  problems.  It  has  proven  to  produce  good 


answers  with  fewer  replications  than  previous  solution 
techniques. 

4.  Successive  Maximization.  It  not  only  maximizes  the  minimum 

distance,  it  successively  maximizes  the  second  minimum,  third 
minimum,  and  so  on.  Because  all  points  and  their  distances  in 
relation  to  each  other  are  important  in  the  spacing  problem, 
this  is  a  more  appropriate  objective. 

Pi sadvantages 

Here  is  a  list  of  the  method's  major  disadvantages: 

1.  Optimality.  The  method  does  not  guarantee  a  globally  optimal 
solution,  nor  does  it  guarantee  even  a  locally  optimal  solu¬ 
tion.  It  just  guarantees  a  "good"  answer  and  it  may  take 
several  runs  to  get  that  "good"  answer. 


2.  Size  Limitations.  The  program  has  been  proven  on  a  problem  of 
at  most  23  points  and  seven  faces.  There  is  no  guarantee  that 
it  will  work  on  very  large  problems.  Computer  time  is  the  key 
resource  here  because  computer  time  required  increases  exponen¬ 
tially  as  the  number  of  points  increase.  It's  up  to  the  user 
to  weigh  the  value  of  his  or  her  computer  time  to  the  value  of 
a  solution. 


VII.  Summary 


I  hope  you've  enjoyed  this  little  trip  into  the  world  of  poly¬ 
hedrons,  colors,  and  the  like.  If  you've  already  forgotten  what  we've 
learned,  here's  a  quick  refresher.  After  being  introduced  to  the  spacing 
problem,  we  learned  of  some  examples  in  the  real  world,  most  notably  the 
color  spacing  problem.  We  were  educated  on  the  theory  behind  the  color 
spacing  problem  and  we  found  out  how  Carter  and  Carter  proposed  to  solve 
it.  Next,  we  learned  how  the  spacing  problem  is  formulated  mathemati¬ 
cally  and  were  able  to  see  some  solutions  to  simple  problems.  We  also 
found  out  how  the  problem  could  be  solved  using  successive  linear  pro¬ 
gramming.  Next,  we  were  introduced  to  a  new  heuristic  algorithm,  very 
similar  to  Carter  and  Carter's,  that  is  designed  to  solve  the  spacing 
problem.  Finally,  we  found  out  how  well  the  new  algorithm  works.  It  is 
now  time  to  identify  our  conclusions  and  some  suggestions  for  further 
research. 

Conclusions 

The  program  works.  That  was  the  main  objective  of  this  effort  so  in 
that  sense,  I  have  accomplished  what  I  set  out  to  do.  Perhaps  equal  to 
that  was  the  objective  of  providing  AFAMRL  with  reasonable  answers  to  the 
color  spacing  problem  that  are  better  than  educated  guesses.  I  believe 
that  the  algorithm  does  that  as  well.  An  objective  of  my  research  also 
was  to  improve  upon  Carter  and  Carter's  method.  Well,  it  is  difficult  to 
compare  the  two  because  time  constraints  prevented  me  from  doing  a  full 
statistical  analysis.  But,  if  you  consider  that  my  method  runs,  is  easy 
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to  understand  and  use,  and  yields  solutions  that  are  comparable  to  Carter 
and  Carter's,  one  must  consider  it  an  improvement. 

To  be  honest  with  you,  I  originally  thought  that  Carter  and  Carter's 
algorithm  was  a  crud"  first  attempt  generating  suboptimal  solutions  which 
could  easily  be  improved  upon.  It  is  a  credit  to  their  work  that  those 
solutions  are  apparently  not  so  crude.  But  it  must  be  remembered  that 
regardless  of  which  method  is  used,  they  both  require  numerous  replica¬ 
tions  to  come  up  with  an  answer  that  is  only  guaranteed  to  be  "good."  So 
there  is  still  plenty  of  room  for  improvement. 

Suggestions  for  Further  Research 

There  are  two  directions  that  subsequent  research  in  the  area  can 
take.  One  is  to  refine  the  new  algorithm.  The  other,  and  most  exciting 
alternative,  is  to  solve  the  spacing  problem  using  some  form  of  nonlinear 
programming.  The  method  described  in  Section  IV  of  this  report  has  some 
definite  possibilities,  but  the  difficulty  is  that  the  constraints  of  the 
color  region  are  not  linear  and  are  very  difficult  to  derive  mathemati¬ 
cally.  I  could  not  solve  for  them  during  my  work  on  this  thesis.  Carter 
and  Carter  describe  the  region  as  "approximately  a  triangular  prism,"  but 
nowhere  in  the  literature  is  there  a  mathematical  formulation  of  the 
boundaries  (2:2937).  This  could  prove  to  be  a  major  stumbling  block  for 
anyone  attempting  to  solve  the  color  spacing  problem  using  nonlinear 
programming. 

Alternatively,  there  are  numerous  improvements  that  could  be  made  to 
my  new  algorithm  as  well.  The  most  worthwhile  appears  to  be  the  use  of 
vectors  and  matrices  to  help  define  possible  new  positions  for  the  two 
closest  points.  The  present  method  disregards  any  alternative  that  takes 


a  point  out  of  the  feasible  region.  A  method  could  be  devised  that  would 
instead  take  the  point  to  the  boundary  of  the  feasible  region  and  still 
consider  it  as  a  viable  alternative.  Indeed,  it  may  be  a  desirable  move 
because  of  the  tendency  of  the  points  to  migrate  to  the  boundaries  and 
corners  of  the  region.  The  trouble  with  this  method  is  that  it  too  only 
works  for  linear  constraints.  Colors  would  have  to  be  spaced  in  the  red, 
green,  and  blue  coordinate  system,  which  has  linear  boundaries,  and  then 
transformed  into  L*u*v*  coordinates  to  calculate  the  distances,  much  in 
the  manner  that  Carter  and  Carter's  algorithm  does. 

Another  interesting  improvement  would  be  to  use  a  concept  called 
"simulated  annealing"  to  help  space  the  points.  In  this  method,  the 
points  would  be  able  to  actually  move  closer  together  in  order  to  even¬ 
tually  achieve  optimal  spacing.  Refer  to  Kirkpatrick,  Gelatt,  and  Vecchi 
(9)  for  more  on  this  possibility. 

Other  less  drastic  refinements  might  be  to  make  the  program  inter¬ 
active,  more  user-friendly,  or  perhaps  include  a  color-naming  subroutine 
that  would  give  the  user  descriptive  names  of  the  colors  selected.  Also, 
improvements  might  be  made  to  the  speed  and  efficiency  of  the  program.  I 
have  thought  of  a  three-phase  system  of  varying  the  maximum  and  minimum 
step-sizes  and  the  step-size  reducing  scheme  to  hopefully  get  consis¬ 
tently  better  answers  more  quickly.  It  would  consist  of  starting  out 
with  large  maximum  and  minimum  step-sizes  and  halving  the  step-size  at 
each  iteration.  This  would  place  the  points  in  a  roughly  spaced  con¬ 
figuration.  The  next  phase  would  have  a  much  smaller  maximum  step-size, 
reducing  by  four-fifths  to  the  same  minimum  step-size.  This  would  refine 
the  points  to  what  one  would  hope  to  be  a  rough  global  optimum  location 
scheme.  Phase  3  would  concentrate  on  further  refining  the  locations  and 


the  accuracy  with  a  very  small  maximum  step-size  and  a  miniscule  minimum 
step-size.  The  composite  result  would  hopefully  be  impressive. 

Yet  another  improvement  would  be  to  selectively  place  the  points  in 
the  corners  of  the  region  then  place  any  remaining  points  randomly  on  the 
faces,  rather  than  randomly  placing  all  points  in  the  interior  of  the 
region  as  an  initial  location  scheme.  This  procedure  would  conform  more 
closely  to  the  observations  made  in  Section  IV  regarding  simple  problems. 

The  biggest  deficiency  of  this  thesis  effort  was  that  the  algorithm 
was  never  fully  tested  statistically  to  properly  compare  its  results  with 
those  of  Carter  and  Carter.  It  is  highly  recommended  that  any  future 
effort  include  extensive  tests  of  this  nature.  I  have  found  the  topic  to 
be  interesting,  important,  and  satisfying,  and  am  sure  that  anyone  else 
who  pursues  research  in  the  area  will  gain  similar  satisfaction. 


49 


Appendix  A.  Color  Transformations 


This  appendix  briefly  describes  the  transformations  from  red,  green, 
and  blue  color  coordinates  to  the  CIE  L*u*v*  color  space.  If  we  let  (YR, 
yG*  yb)  be  the  luminances  for  the  CRT's  red,  green,  and  blue  guns, 
respectively,  then  the  first  task  is  to  find  something  called  the  tri- 
stimulus  values  (Xy,  Yy,  Zy).  This  transformation  is  given  by  the  fol¬ 
lowing  matrix  operation: 


V*R  *6^6  V^B 
1  1  1 

zr/-Yr  zG/yG  ZB^B 


where 

(xR,  yR,  zR)  =  the  chromaticity  coordinates  of  the  red  gun 

(xq,  yg,  Zg)  =  the  chromaticity  coordinates  of  the  green  gun 

(xg,  yg,  zB)  =  the  chromaticity  coordinates  of  the  blue  gun 


These  chromaticity  coordinates  correspond  to  the  locations  on  a  chroma¬ 
ticity  diagram  of  the  red,  green,  and  blue  colors  corresponding  to  that 
particular  cathode-ray  tube  gun,  and  are  values  between  zero  and  one.  In 
addition,  x  +  y  +  z  =  1  for  each  gun. 

The  next  step  is  to  find  the  u '  and  v'  values  for  each  color.  These 
values  correspond  to  the  color's  location  on  the  CIE-UCS  chromaticity 
diagram  shown  in  Figure  2.  The  transformations  are  given  by: 


"V  'T  ■  l  >  u  •»  ■  '  '  ■  "AJ 


U>  I^iyij  Iiiitfii  ■  «.•  J  u.«  t,"  ‘  -  ■  '  V  *T 


u'  =  4x/(-2x  +  12y  +  3) 

v"  =  9y/(-2x  +  12y  +  3) 


where 

X  =  Xj/ ( X  y  +  Yy  +  Zy ) 
y  =  YT/(Xy  +  Yy  +  Zy) 

The  final  step  is  to  convert  from  (Yy,  u',  v")  to  (L*,  u*,  v*)  coor¬ 
dinates.  This  is  accomplished  by  the  following  set  of  equations: 

L*  =  116  ( Y y/Y q ) 1 /3  -  16 

u*  =  13  L*  (u'  -  Ug) 

v*  =  13  L*  (v'  -  Vg) 


where 

Yg  =  the  value  of  Yy  when  all  three  guns  are  turned  on  full 
blast,  producing  a  pure  white  color 
Ug  =  the  value  of  o'  for  pure  white 
Vg  =  the  value  of  v'  for  pure  white 

That’s  all.  Pretty  simple  isn’t  it?  The  reverse  transformations  are 
equally  as  simple  and  are  contained  in  Appendix  3  (which  is  just  a 
listing  of  the  program)  among  the  coding  of  subroutine  "RG80UT." 
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Appendix  B. 


FORTRAN  Code  and  Doe unn-n tat  ion 


( ;  #  *  X  »:  ;«  *:  t  >j;  *  A  x  *  t  y/  >!■:  >K  *  #  $  %  t  .A ;{;  :  j(c  t  >|e ; 


»A-  •«/  \V  \V  *.'  ■  viy  \V  V/  Vj-  'A  vV  V/  \V 

'l-  l-  ,  '(  'A  'J-  •  '(-'I-  ’I-  'i-  'I-  'I • 


(.: 

(.:  nil  ms  program 

l  BY 

C  ROBB  cm  EY 

C  00R--B4D 


c 

C;  *  ^  >K  >K )« >J.  >{<  >!-: 


He  %  X  X  t He  *  X-  X  He  He  He  He  He  He  H<  H<  He  He  He  He  He  *  He  He  He  He  He  He  He  He  He  He  He  He  He 


He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  He  H  He 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c: 

c 

r: 

c. 


HA IN  PROGRAM 

THE  PUR ROBE  OF  THIS  PROGRAM  I B  TO  MAXIMIZE  THE  MINIMUM 
D I  STANCE  AMONG  N  POINTS  IN  A  CONVEX  THREE  -111  ME  NS  I  ON  AL  REGION 
AND  THEN  SUCCESSIVELY  MAXIMIZE  THE  SECOND  SMALLEST  DISTANCE, 

THE  THIRD  SMALLEST,  AND  SO  ON,  UNTIL  AL.L  POINTS  ARE  OPTIMALLY 
SPACED*  THE  REQUIRED  INPUTS  FOR  THIS  PROGRAM  ARE  THE  NUMBER 
OF  POINTS  BEING  SPACED  <N),  THE  DESIRED  INITIAL  STEP-SIZE  (SIZE) 
THE  NUMBER  OF  CONSTRAINTS  (FACES)  ON  THE  REGION  CF),  AND  THE 
EQUATIONS  OF  THOSE  CONSTRAINTS  IN  THE  FORM  AX  +  BY  +  CZ  <  K  . 
THE  PROGRAM  THEN  OUTPUTS  THE  OPTIMAL  VALUES  OF  X,  Y,  AND  Z, 

AND  THE  FIRST  TEN  (OR  SOME  OTHER  NUMBER  CHOSEN  BY  THE  USER) 
MINIMUM  DISTANCES* 


NOTE?  THROUGHOUT  THIS  PROGRAM,  X,  Y,  AND  Z  ARE  USED  TO  DENOTE 
THE  LOCATIONS  OF  THE  POINTS  BEING  SPACED  AND  CORRESPOND  TO 
LJK ,  U#f  AND  V#  FOR  THE  COLOR  SPACING  PROBLEM.  THIS  SHOULD  NOT 
BE  CONFUSED  WITH  THE  TRISTIMULUS  VALUES  OF  A  COLOR  WHICH  ARE 
COMMONLY  REFERRED  TO  AS  CAPITAL.  X,  Y,  AND  Z  IN  THE  TRADE.  IN 
THIS  PROGRAM,  THE  TRISTIMULUS  VALUES  WILL  BE  REFERRED  TO  AS 
X TOTAL,  Y TOTAL,  AND  7. TOTAL* 


THERE  ARE  SEVERAL  NOTES  TO  THE  USER  IN  THE  PROGRAM  WHICH 
ARE  USED  TO  SIGNIFY  PLACES  WHERE  YOU  MAY  BE  REQUIRED  TO  MAKE 
CHANGES.  THEY  ARE  SUMMARIZED  FIERI::  FOR  YOUR  CONVENIENCE* 


IN  THE  MAIN  PROGRAM,  THE  USER  IS  RF  OUT  RED  TO  SUPPLY  Jill 
NUMBER  OF  POINTS  TO  BE  SPACED,  THE  SLED,  IMF  MAXIMUM  STEP-BIZ 
THE  MINIMUM  STEP-SIZE,  AND  TUT  STEP-STZE  REDUCING  SCHEME.  IN 
SUBROUTINE  PEACE  POINTS,  HE  MUST  SUPPLY  HIE  MAXIMUM  AND  MIN¬ 
IMUM  VALUES  ALLOWABLE  FOR  EACH  01  Till:  THREE  COORDINATES.  IN 
SUBROUTINE  FIND  B  C STANCE ,  Ti  ll  USER  MUST  CHOOSE  THE  NUMBER  OF 
DISTANCES  HE  WISHES  TO  APPEAR  ON  Jill  PROGRAM'S  OUTPUT. 


, 


C 


I  N  S I J  «  R  0  0  (INF  ('ll  I.  0 1\  -  CO  I  N  I  ,  I  HI.  0  S I  I  \  I  S  R  F  0  0  T  R I  1 1  l  0  S  LJ  P  PLY 


r)2 


i  :  ,  ■ .  ■  i :  ,i  ■  : : ‘  1 ;  1  I  i  •.  1  ■ 1 .  A •  i i  ■  l  H>  t  i > . 1  A  i  i  (if) 

c  |  i  ■ :  ill!.'.:  !  !!!!  I  In;.,  MX  i  IP  F  !.  /  :  I. ,  I  01:  10! 

I'  COLOR  '  ! •  A  |  :  l-KOfU  I  ii,  Nil.  NUMBER  0!  FACES  SIIOUI.  0  USUALLY 

(j  I-t|  KEVIN.  IN  A  Tip  t  T  l ! )  N  i  All.  0  II II  1%  VARIABI  LS  STAY  TUT  SAMI" 

I  EXCEPT  r  01 '  111!  VAI  UPS  OF  K;  K(1)  is  thf.  maximum  aelowabi  f 

('  I  UfilNAMLI.  VO"  l.lp.ll  I"  SI  I  I  OL  nil  Kin  (•OK1  v  K  ( 2 )  -  THF.  MAXIMUM 

L  Ills  EN  1.1  if-!  I  U  ■  i  i A f > ;  !  ..!•-.)  1 1 II  MAXIMUM  III.  OF  01  IN  IUHTNANCE, 

0  AMII  !•  '  '  •  I  HI  ill  I'M  mom  !  O  I  Al  lUMINANOI  YOU  WISH  10  PF  KM  I  1  (IT 

i;  s 1 1 o o I  1 1  pi  i  i-m  I,  A"  a  on. a i  i or:  value). 

C; 

r  |  0  Si  i ' : ! : !  r  I  •  :  IF!  I  |.  I'M  I  0  t  III!  ONI  V  III  "( I  I  III  0  CHAKOL  IB  T  III 

o  val  or  or-  i  n;  rj, r  mom  s  i  i  ;  •  six,  • 

c 

C  1  O':  (  III  UK  SPACING  PROBI  KHSV  BOBKO 01  INF  ROIUJOT  IS  VERY 

0  [Ml  'OR  FAN  I  <-  II I  FRF  ,  IT  IF  USER  MUST  SUPPLY  IMF  Cl  IK  DM  AT  [CITY 

n  COORD  IN AT  I:  S  OF  1 1 II..  BUMS,  THF  MAXIM IM  I.UMI  KAN  Cl  ALLOWABLE: 

C  (I -HOT) ,  AND  THE  VALUES  or  U  PRIME  MOT  AMD  0- PRIME  MOT,  THAT 

C  IS  Till.  I  XT  I  Ml  OF  I  HE  CHANGES  I  HAT  NEED  TO  BE  MADE  I- OR  ANY 

C  NORMA!  USE  OF  I  II  [S  PROGRAM. 

C 

C 


C 

C  OAR  l  ABLE  DE I  I N I T I ONS 

C 

C  real,  variables: 

C  X(I)=  A  0  N  E  - 1. 1 I M  P.  N  S 1 0  N  A  l...  ARRAY  CONTAINING  THF  /.-COORDINATE  OF 

C  POINT  I 

C  Y(I>~  THE  Y-  COORDINATE  OF  POINT  I 

0  Z  < I > ~  THE  7-COORD T MATE  OF  POINT  I 

(I  D(  I  ,,))-•  THE  EUCLIDEAN  DISTANCE  BETWEEN  POINT  I  AND  POINT  J 

C  CAPD-  THE  CURRENT  MINIMUM  D  [STANCE  BETWEEN  ALL.  PAIRS  OF  POINTS 

C  NEWD-  THE  NEW  MINIMUM  DISTANCE 

C  SIZE-  THE  STEP-SIZE 

C  IFNDX  ( I )  -  THE  X  -COORD  I  NATL  OF  ENDPOINT  I,  WHERE  THE  ENDPOINTS 

C  ARE  THE  TWO  CLOSEST  POINTS  IN  THE  REGION 

C  FNDY < I } "  THE  Y  COORDINATE  OF  ENDPOINT  I 

C  E  N  D  7.  <  I'  )  -  1 1 1 F  Z  •  -  C  0  0 1  s  0 1 M  A  I  E  0  F  E  M  D  P  0 1 N  T  I 

C  All  X  ( I ,  J  )  -  (TIE  ,  IT  IT  A I  TERN  A  IE.  I..OCAT.TON  FOR  THE  X~  COORDINATE  OF 

C  ENDPOINT  I 

C  AL.TY  ( I ,  ,1 )  ---  THE  JTH  AL  TERNATE  LOCATION  FOR  THE  Y-COOROI NATE  OF 

f  ENDPOINT  I 

C  AI..TZ  (  I ,  J ) :  THF  JTTI  A I  TERN  ATE  LOCATION  FOR  FITE  7-COORDINATE  OF 

C  IN  lip  (TIN  I  I 

C  Ml  .MX  (  I  >  1 1 1 1  X  ••  COORDINATE  OF  POINT  I  PROM  I  HE  PREVIOUS  ITERATION 

C  ME  MY  v  I  )  THE  \  COORDINATE  01  10  [.NT  I.  FROM  THE  PREVIOUS  ITERATION 

C  Ml.  HZ  (  I  )  THF  Z  COORDINAT  E  0!  POINT  I  FROM  THF  PREVIOUS  ITERATION 

C  DUMJ, DUMP, DUMB-  DUMMY  VARIABLES  REQUIRED  FOR  SUBROUTINE  RGBOUT 

C  BUT  NOT  USED  IN  THE  MAIN  PROGRAM 

C  SPED  HIE  STARTING  VALUE  FOR  THE  RANDOM  NUMBER  GENERATING  SEQUENCE 

C  IN  THE  CYBER  KNOWN  AS  SUBROUTINE  ' RAN SET ' 
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I 


N  ! :  ■;  ;i  ii'ii  :  I ;  n;  i  ;  i  | : :  i  1 .  i i  >i  i  i  : ;  i  m .  .  .1  Y, i  |  tH> 

I  NDi  I  I  III!  (.'Uril'-M;  AS/-1  H.  !  n  1 1  P  l-l  III!  Ii."iROINI  ! 

I  I  i  i !  i ; '  I  '  -III  i-llMlii  !:  A'  SHI  IAI!  u  W  I  III  !  N 1 U  '()  I  Ml  2 

C 

i  l  i  iai  :••,[  ii  i :  uai  •  I  ri Hi.  I  s  ; 

<  S  T  A 1 1 ' '  >  (III  AIAH.IS  01  111:  RROURAM;  'RUN'  II  I !  !t.  STOPPING. 

i  f  ii :  i:  rni-)  has  mu  i.!  pa  mi,,  and  •  done  •  if  it  has 

I  ' .  I  r'i  i  ; '  I  i  '  r  I  :i  SIAIUS  (i|  I'M!  j\i  1  i  ;  '|  |(IR  |  |  XI.H  v  AND  'IF' 

I  :  Ml  .  !  ,;■!!  I  '  i  II 

r 

r.  I  Mil  REAL  ‘IAI-:  I  ADI  i  St. 

C  FAN  I  '.'Al  Si  ADI.  I  HSi  li  IN  I  III  FANDOM  I'UMCIF  FI'  f  1  If:  RAT ’.I  KG  SUB  ROUT  I  NF 

F  IN  I  III  F  r'DI  I:  (.  Al  I  I  D  FANS!  i 


i: 

(  ;y  ;V  y.  >y  ;v  j-  ;.i.  |  g,  g,  ■)/  g/  g  g.  g,  1/  vi,  g 
F 


t  %  1 Z  Z  I' Z :{:  Z.  t.  %  >;  t  >!. 


*  Z  Z  Z  Z  Z  >:■;  Z  Z  t  Z  #  *  He  %  Z Z 1 Z  Z  Z  t  %  Z  #  Z-  H 


C 


PROGRAM  SPACP ) 

R  F  A I  X  <  50 )  ,Y<  5  0  )  v  Z  <  5  0 )  >•  D  < 5 0 .5  0  )  ,■  C  A  P D ,  N  E  W  D  ,  B I Z  F  ,  E  N  D  X  <  2 )  ,  E  NDY<2)  , 

$  E  N  D  Z  ( 2  )  » A  FIX  (  2 , 2  7  )  ,  A  L  T  Y  (  2 ,27)  ,A  I...  T  Z  <  2 ,27),  M  E  M  X  <  5  0 )  ,  M  E  MY  (BO) 

X  M  L-.  M  Z  (50),  D  U  M 1  (  5  0  )  ,  I  UJ  M  2  <  5  0  )  ,  D  U  M  3  (  5  0  )  ,  B  E  E  D 

I N T E G E R  N  ,  END P T  I.  ,  EM D P T 2 


C I  I A R  A C  T  E  R  ST  A  T  U S  Z-  4  ,  B  T  A  T  P  1(5 0 )  *  I 

EXTERNAL  RANF 

OPEN  (  3  ,  F I L  E~  SPG  OUT  '  ) 

REWIND  3 


10  "3 
I!  DO  "3 


c  . -  IMPORTANT  NOTE  TO  USER  -- 

C 


C  THE  USER  MUST  SUPPLY  THE  VALUE  OF  N  AND  SIZE  HERE,  AND  THE  SIZE 

C  AGAIN  IN  SUBROUTINE  '  F IX-  POINT , '  I  RECOMMEND  USING  A  STEP-SI  ZE 

C  APPROXIMATELY  EQUAL  TO  ONE -HALF  THE  DISTANCE  BETWEEN  THE  TWO 

C  FARTHEST  POINTS  IN  THE  REGION,  BUT  THE  USER  IS  ENCOURAGED  TO  TRY 

0  THE  SAME  PROP!.. EM  WITH  SEVERAL  DIFFERENT  MAXIMUM  STEP-SIZES  TO 

C  INCREASE  THE  CHANCES  OF  GETTING  AN  OPTIMAL.  SOLUTION. 


SIZE- 130.0 
N~  10 

S  T A  i  ns;  •  '  DONE  ' 

DO  1.00  !  I.  v  N 

STATPTC  f )  -  ' U ' 
it  I.TT  X  (  T  >  0 . 0 
ME  MY  (  .[  )-().() 
MLMZ  (  I.  )  -  0.0 
1.00  CONTINUE 


IMPORTANT  NOTE  TO  USER  - . 


o  n  n 


I  1:1  .  ..  1  ;  "  '  .  ;  -  .1  :  i  i  .  i .  Vf.i.  Hi  |(i),  |||;  LANDHC  PUn|;l  !'  I  '  I  '■ .  i  i  ! ■ 

l  81..  I H  It  Ml  i  iailJi.il  RANDOMLY  CIIOOSLS  I  III  INI. THAI  LOLA! IONS  181,  I  I  i: 

c  i'(irN!8.  :n  1 8  in  commlnui  d  nir.i  h  ii  n  lit.  rums  Bi  pi  :m  okmi  d  i  or 

i  ai  i  pRiii'i.  i  h ■ > ,  in  miii  ci i  carp  mi  user  8II0UI  n  i: hanoi:.:  tip  sit  u  for 

c  1  A  Cl  I  I-:  (.IN  ,  IT  CAM  BI  ANY  N  I  N!  HI  (PIT  IH1LGLR, 

C 

81  I- 1 1  .'I'.VCHA!  21  ,0 
UR  I  TF  •:  A  ,  J.  I  0)81,  n 

I  1 0  r  0 1: M A T  (  I  6 X  v  'Till  811  D  : '  ,1  1  I  -  I  > 

i:  THIS  RFC  I  I.  ON  CAI..(  0:  n  i  PR  TUP  INI  T  in  I  LUCA  I  I  Oil  8  I  OR  :  III  POINTS  AND 

C  F  0  R  H  A  I :  A  N  T  i  P I  •:  L  N  T  8  III  T.  I...  0  C  A  1  I  0  N  8  . 

CALL  RANOET(BFFD) 

C  A I ...  I...  P I  A  C  P  T  <  N  y  X  Y  v  Z  r  8  E  Li  D ) 

MR I T  E ( 3 , 1 2 0 ) 

J  20  FORMAT ( 11 X THESE  ARE  THE  STARTING  LOCATIONS') 

CALL  PRINT ( X  ,  Y  ,  Z  ,  N  ) 

C  A  I. .  L.  R  G  B  0  U  T  (  X  ,  Y  ,  Z  ,  N  ,  D  U  M  .1.  ,  D  U  M  2  ,  D I J  M  ,5 ,  S  T  A  T  U  S  ) 

STATUS- ' RUM ' 

C  THE  MEAT  OF  THE  PROGRAM  BEGINS  HERE 

1 3  0  CAL.  L  FINDS  <  X  ,  Y  ,  Z  ,  N ,  S  T  AT  FT  ,  C  A  P  D  ,  E  N  D  X  ,  E  N  D  Y  ,  E  N  D  Z  ,  E  N  D  P  T 1 ,  E  N  D  P  T  2  , 

8  STATUS) 

C  A  L  L  D  Ii!  F  ALT  ( E  N  D  X  ,  E  N  D  Y  ,  E  N  D  Z. ,  S I Z  E  ,  E  N  D  P  T 1 ,  E  N  D  P  T  2  >•  S  T  A  T  P  I  , 

X  ALTX , ALTY , AL TZ ) 

CALL  CHKALT  <  X  ,  Y  ,  Z  ,  N CARD  ,  ALTX  ,  ALTY  ,  ALTZ  ,  STATPT  ,  ENDPT 1 ,  ENDPT2  ,  NEWD  ) 


C  IF  THE  MINIMUM  DISTANCE  IS  GREATER ,  THEN  KEEP  ITERATING 
IF  (NEWD  .  GT.  CARD)  THEN 
C A PD “NEWD 
GO  TO  130 

C  IF  NOT,  THEM  REDUCE  THE  STEP -SIZE 
ELSE 

SIZE.” SIZE/ 2 . 0 
IF  (SIZE  >LE>  ,0001)  THEN 

NOTE  J  THE  USER  MAY  WISH  TO  USE  A  DIFFERENT  MINIMUM  STEP- SIZE, 

FOR  THREE  DECIMAL  PLACES  OF  ACCURACY,  I  USE  ,0001.  YOU 
MAY  AI..SO  WISH  TO  USE  A  SCHEME  OTHER  THAN  HALVING  THE  STEP- 
C  SIZE,  LIKE  REDUCING  IT  BY  THREE -FOURTHS  INSTEAD. 

CALL.  FIXPT  (  X  ,  Y  ,  Z ,  N ,  STATPT  ,  ENDPT  1 ,  ENDPT 2 ,  SIZE  ,  STATUS  , 
a  MEMX , MEMY , MEMZ ) 

I F  (  S  r  A  T  U  S  ,  E  Q  ,  '  R  U  N '  )  T  H  E  N 
GO  TO  130 

E L  S E  I F  ( S T  A T U S  . E 0 ,  ' D 0 N E ' )  THE N 
GO  TO  1  TO 
END  IF 

ELSE 

00  TO  130 
END  IF 


END  IF 

C  THIS  SECTION  CONTROLS  THE  FORMAT  OF  THE  RESULTS 


L  TO 
150 


WRITE (3, 150) 

FORMAT ( 15X, 'THIS  IS  THE  FINAL  ANSWER') 

CALL  PR I NT (X,Y ,Z,N) 

C A L. L  R G BOUT  (  X  ,  Y  ,  Z  ,  N  ,  D U M .1.  ,  D U M 2  ,.  D U M 3  ,  S T  A T U S  ) 


-AcVA1* 

Cot* .  V,.YH  ‘ 


\o 


n 


,vc 


&o ' 


yo\ 

; 
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c 

c. 


7 1  ii  I  'URI  '  IS!::  I  ll 

in  tun  r<r  n t on,  the 


1 1!  IS  BABY  IS  TO  RAN  ITCH  I  Y  I'l  A  (II:  Till.  POINTS 
LITTLE  TYKE  RECEIVES  THE  NUMBER  OF  POINTS 


( N )  A  8  AN  INPUT  ,  AN  IT  GURGLES 


1 1 1 1  B  T  A  R  ?  .1  N  G  I  0  C  A  I  .1  0  N  B  !:  0  R  Till 


POINTS  (X,  Y ,  A  TED  X)  AS  OUTPUTS. 


C  NOTE  ?  THIS  BUB  ROUT  INE  OBI'S  THE  RANDOM  NUMBER  GENERATING  FUNCTION 

C  IN  THE  CYBER  CALI. ED  '  RANSET .  '  MOST  MAIN -I- EAME  COMPUTERS  HAVE 

C  A  RANDOM  NUMBER  GENERATING  CAPABILITY,  BO  CONSULT  YOUR  LOCAL 

C  COMPUTER  EXPERT  FOR  DIRECTIONS  ON  HOW  TO  USE  IT  ON  YOUR  FAVORITE 

C  SYSTEM. 

C 

C 

C 

C  .  VARIABLE  DEFINITIONS 

C 

C  REAL  VARIABLES*. 

C  X(I)=  THE  INITIAL  X -COORD INATE  OF  POINT  I 

C  Y  C I  > -■=  THE  INITIAL  Y -COORDINATE  OF  POINT  I 

c  z  <  i  )  =  t i  i e  :r. n i t i  a l  z - c o o r d i n a t e  of  p o i n t  i 

C  YRED  ( I  ) THE  RED  PHOSPHOR  LUMINANCE  FOR  POINT  I 

C  Y GREEN <  I ) ::::  THE  GREEN  PHOSPHOR  LUMINANCE  FOR  POINT  I 

C  Y  B  I..  IJ  E  ( I )  =  T  H  E  B  I..U  E  P  H  0  S  P  H  0  R  L  U  M I N  A  N  C  E  F  0  R  P  0 1 N  T  I 

C  SEED--  THE  STARTING  VALUE  FOR  THE  RANDOM  NUMBER  GENERATING  SEQUENCE 

C  IN  THE  CYBER  KNOWN  AS  SUBROUTINE  'RANSET' 

C 

C  INTEGER  VARIABLES: 

C  N~  THE  NUMBER  OF  POINTS  BEING  LOCATED 

C  J~  THE  NUMBER  OF  POINTS  BEING  TRANSFORMED  BY  SUBROUTINE  '  RGBOUT ' 

C 

!.  CHARACTER  VARIABLES : 

C  POINT-  THE  LOCATION  OF  THE  POINT,*  EITHER  'IN'  OR  'OUT'  OF  THE 

0  REGION 

C  DUM.E  A  DUMMY  VARIABLE  REQUIRED  FOR  SUBROUTINE  'RGBOUT'  BUT  NOT 

c  ubf:d  in  this  subroutine 

(.: 

C  I  XTI  RNAI.  variadi.es: 

c  RANI-  A  VARIABLE  USED  IN  CONJUNCTION  WITH  'RANSET  '  TO  GEN!  RAT  I 

C  A  RANDOM  NUMBER  BETWEEN  ZERO  AND  ONE 


C 

c 


**  *:#  *  %  #  *  %  *  *  *  *  **  *  %  *  %  %  %  %  %%  %  %  %  %  %%  %%  %  %  %%%  *  %%  %  *  %  %  %%  %  %  %%  Hoy  X  #  *  *  #  *  y* 
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I  .'0 


!  1  mi.  i  ■  i  i  .i  r  \  <  m  ,  x  ,■  •.  -  v ,  m  i  u  > 

I'! h:  A I.  X  ( S 0  i  v  Y  ( ' 0  >  ,  /  ( I.  i  ll  v  Y I  :l  II  ( 0 0  >  Y 0 R I.  I .. N  ( 0  )  ,  Y I H . 0 L  ( 0 < > ) 

I.NTFGI  I.  N  ,  J 

C I  I A I  ■:  A C  (  1. 1 :  I '  0  l  j  *  I X  ,  T 1 1 1 H  |  *  4 

I  X  T  FRwAI  RANI 
DO  1  AN  l  l  ,  N 

X  (  T.  >  :  RANI  (  ) (  I  ON  .  0-7 . 0  )  {  V  ,  0 

y  (  i  >  i •: a n i  (  >  x  ( i ?x .  o  i  .  o  > .  o 

x  (  i  >  r  a  ;■  i  r  c  )  *  (  i.  o  o  ,  o  i  i  7 s ,  o'  >  1. 2 s .  0 


C  --  I  Ml  'OR  TAN  I  NO  I 1  Id  URL  R  . 

r; 

C  I  I II:  VARIIARI  L  •RANI''  GIVES  VAl  dl  R  IN  I  Hi  RANGE  FROM  ZERO  TO  ONI:  . 

C  IF  IMF  USER  WANTS  VALUES  0(1 1 S  T  HI  OF  THAT  RANGE,  THE  EQUATIONS 

C  FOR  X  y  Y ,  ANTI  Zy  SHOULD  I  OOF  I  IK  I  THIS: 

C  X  (  T  ) ~  R  A  N  F  (  )  *  (  X  M  A  X  X  M I. N  )  T  X  M I N 

c  y ( :r. ) - r a f i f'o * •: v m a x y m :i n )  +• y m i n 

t z  <  :i:  >  -  r  a  n f  ■;  >  *  <  z  m  a  x  -  z  h  i  n  >  +  z  m  i  n 

C  WHERE  MAX  AND  MIN  ARE  IT  IF  MAXIMUM  AND  MINIMUM  POSSIBLE  VALUES 

C  FOR  X,  Y,  AND  Z,  RESPECTIVELY. 

C 


CALL  RGBOUT  <  X  ( I.  >  ,  Y  <  I  )  ,  2  ( I  >  ,  J ,  YRED  <  I )  ,  Y  GREEN  (I )  ,  YBLUE  (I )  ,  DOM  I  ) 
C  A  I...  I...  C  H  E  KPT  (  Y  R  ED  <  I  >  ,  Y  GREEN  <  [  >  ,  Y  B  L  U  E  ( I )  ,  P  DINT) 

IF  (POINT  .EQ.  'OUT')  THEN 
GO  TO  170 
END  IF 

1.60  CONTINUE 
END 
C 


C 

C  %  y.  >!<  *  %  y.  %  y  >!•:  %  %  %  %  %  %  %  %  %  %  :{<  %  %  %  %  *  #  #  1 1 1  # ;fc  *  #  1 t y.  *  y.  #  :)< 

C 

C  S U  B  R  OUT  I N E  F I N B - D I S TAN  C E 


C 


C  THE  PURPOSE  OF  THIS  SUBROUTINE  IS  TO  FIND  THE  MINIMUM  BIS- 

C  TANCE  BETWEEN  AU...  N(N~:L  )/2  PAIRS  OF  (UNFIXED)  POINTS.  THE 

C  INPUTS  FOR  THIS  PROGRAM  ARE  THE  LOCATIONS  OF  THE  POINTS  (X, 

C  Y,  AMD  Z),  THE  STATUS  OF  THOSE  POINTS  AS  TO  WHETHER  THEY  ARE 

C  FIXED  OR  UNFIXED  (STATPT),  AND  THE  NUMBER  OF  POINTS  (N),  THE 

C  PROGRAM  THEN  OUTPUTS  THE  MINIMUM  EUCLIDEAN  DISTANCE  BETWEEN 

C  ALL  PAIRS  OF  (UNFIXED)  POINTS,  THE  TWO  POINTS  THAT  CORRESPOND 

r  TO  THAT  MINIMUM  DISTANCE  ( ENDPT  I  AMD  F.NDPT2)  ,  AND  THE 

0  COORDINATES  OF  THOSE  POINTS  (ENDX,  I  INDY y  AND  ENDZ) . 

(.:. 

c 

C 

L'  V  A  R I  A  B  1.1.  D  1:1 1  I N  I  T  1 0  N  S 

.  c: 

c  RI.IAI  variables: 
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I 


I 


c 

/  •:  !  i 

1 ;  i :  .■■■■  1  1 1 1 ;  i 

TINA  II  in 

Id!  N  I  1 

CAPO 

i|  IF  MIN  IMS 

iii  HU  I  1  HI  1, 

N  IT  1  A  I  Al  M  I. 

BL.TIJI  I  N  Al 

1  ( IJN!  i  XI.  D  )  P 

c 

HI  ITU  NTS 

L 

1.  NT  1-  (  1 

>  III' 

■(TURD  1.  NAT  1 

Ill  |  j-1  (Till  1  i ! 

IT  1 

f' 

UNITY  <  1 

■  nn  v  f 

l  TIED  1  N A  1  L 

01  .  r.THTi  1  f)  1  1 

C 

L  ND'/  (  1 

1  hit  /  i 

THIRD  1 1-) A  1  1 

III  1  -'Hi  Till 

1  I  1 

c: 

!T(  1  ,  ,  1  > 

1111  in  c:  1 

IIU  AN  HI  Si 

ANi'L  1  fill  11 

I'O  INI  1  10 

ITT  1  NT  ) 

G 

\" 

it  1  1  i  . 

1  III  !  HI 

SMnl  .i  1  ' .  i  D 

l  S  r  i  ,■  H  ; 

1 . 

c: 

1  NT  LOP 

R  Uni'Ll  ABLE 

».>  < 

c 

N  I  I  IF 

NUMBER  OL 

POINTS 

c 

ENDPT! 

-  Till:  NUMBER  ASSOC  I A 

I  ED  Will-!  E 

N DP 01 NT  1 

n 

ENDPT  2 

-  THE  NUMB 

1  R  ASSOCIATED  WITH  E 

NDPOINT  2 

c 

K -  T  I  II 

COUNTER  U 

IS  ED  IN  ORDERING  Till: 

DISTANCES 

c 

L  THE 

NUMBER  OF 

DISTANCES 

THE  USER 

WISHES  TO  T 

•'PINT 

c 

EYE <  K ) 

-  AN  ENDPOINT  CHORES 

PONDING  TO 

THE  Kill  Si 

TALLEST  DCIvJ) 

c 

JAY ( K ) 

»  THE  OT HE 

E  ENDPOINT 

COORESPONDING  TO  THE  KTH  SMALLEST 

II  ( .1: ,  J ) 


C  C  l-l  A  Is  ACT  E  E  K>  A  IS  T  A  B I ..  E  S  *. 

C  STATPT  CD-  THE  STATUE  OF  POINT  I?  EITHER  FIXED  OR  UNFIXED 

C  STATUS-  THE  STATUS  OF  THE  PROGRAM;  EITHER  'DONE'  OR  'RUN' 

C 

c 

c#******#*#*  :i<  *  #  *  »:  *  *  #  lI<  *  $  :{<  *  %  t  *  >}'  *  >K  **#  *  ***  %  *  sic#**#*#*######  X<  ####  *  *  *  *  *  *  *  *»<  #  * 

C 

C 

S  U  D  R  0  U  TIN  E  F I N  D  D  ( X  „  Y  ,  7  ,  N S  T  A  T  P  T  C  A  P  D  ,  E  N  D  X  ,  E  N  D  Y  ,  E  N  D  7  r  E  N  D  r:'  T 1 P  E  N  D  P  T  2  v 
&  STATUS) 

R  E  A  L  X  <  5  0  )  ,  Y  (  5  0  )  YZ<  5  0 )  C  A  P  D  ,  E  N  D  X  (  2  )  E  N  D  Y  <  2 )  ,  E  N  D  Z  <  2 )  ,  D  (  5  0  r  5  0  )  „ 

8  D IS (300) 

I N TEG E R  N  ,  E N D P T 1 E N D P T 2 K ,  L  ,  E  Y E  (  3 0 0  )  ,  J A  Y  (  3 0 0  ) 

CHARACTER  STATPT  ( 50 )  *  j. ,  ST  AT  US*  A 
GAT'D”  1 0000  <■  0 
DO  200  I  ~  .1.  ?  N  •  J. 

DO  210  J-I  +  .1  ,  N 

C  UNLESS  BOTH  POINTS  ARE  FIXED,  CALCULATE  THE  DISTANCE  BETWEEN  THEM 

I F  (  (  STATPT  <  i:  )  ,  LO  ,  '  U  '  )  OR  .  (  ST  A  TPT  (  J  )  .  EC! ,  '  U  •'  )  ) 

8  THEN 

C  THIS  IS  THE  EUCLIDEAN  DISTANCE  FORMULA 

D  ( I  r  J )  ~S  OR  T  C  (  X  C I  )  X  (  J  )  )  **2-E  (  Y  ( I ) -Y  <  J )  )  *  *  2  L 


(  7.  (  I  )  -7.  <  J )  )  #  Y  2  ) 

IF  (  D  C  I  ,  J  )  .1  I  ,  CAI  'Li )  THEN 

CAP  D-D <  I , J) 

END  P  L  l.~l 
e:n  DPT  2  I 
END  IF 


END  IT- 

210  CONTINUE 

200  CONTINUE 

ENDX ( t ) ~X ( ENDPT 1 ) 


CopT 


n 


W 


\o 


7 


fiC 


Ac- 


pc' 


vnP- 


T\0' 


58 


I  .Mil  X  ( 2  1 : "  X  <  I:  N  l.i  r  T 2 ) 
i  N T i Y  <  2  >  V  < !  NDI  ’12) 
I  NDZ  ( •’  >  ;•  < I  i  * ( t i  •  T  2  ) 


nils  si: (  i  :i  i.in  ni : « •  i  i : ' .  1  hi  distances  and  lrin  i s  mi:  i.  smallist 
i.i  i  stan: 'i:  3  in  nu  r;iii  win  n  nit;:  nr:  ns  ram  is  done 
!l  (ST  All 'S  <  I  (S  TK.iNI.  '  )  II  UN 
I... "  1 0 

HI!:.  1 1  Si:  R  CUNT  KOI  .S  II  li  N 1!  i'i  B  1: 1 :  01-  HI  SI  ANCES  HE  WISHES  TO  SI-  II  BY 
CHANGING  THE  VALUE  01  I.  HERE,-  HE  CAN  CHOOSE  TO  SEE  OP  TO  TOO 
H  I  STANCES , 

DO  220  l\  ■  1  v  TOO 

D  I  S  (  K  >  '"1000 . 0 
)  CONTINUE 

DO  230 

DO  2 TO  J~I  f 1 y N 

DO  250  K -  I  ,  I.. 

:i:  f  ( r.i  <  r ,  j  >  .  l  t  .  d  :i.  s  <  k  )  >  i  i-i  e  n 

DO  260  M-L y  K  +  l f -1 

d:i'S(M)-d::s(m-i  ) 

EYE  (  M )  -EYE  <  M- 1 ) 

JAY  (  H ) -JAY  ( H-l ) 

)  CONTINUE 

DIS(K)-IK  I , J) 

EYE  <  K )  ■=  I 
JAY (K ) ~ J 
GO  TO  240 
END  IF 

)  CONTINUE 

)  CONTINUE 

)  CONTINUE 

WRITE <3y280)L 
DO  2 VO  K •“  J.  ,  l. 

Ul  R I T  E  (3,2  V  5 )  E  Y  E  (K  )  ,  ...1 A  Y  (  K  >  ,  D I S  <  K  ) 

)  CONTINUE 

WRITE (3 ,270) 

)  FORMAT ('  ') 

)  FORMAT (VX, 'THESE  ARE  THE  M3,'  SMALLEST  DISTANCES') 

,  FORMAT < 1 VX ,  ' D ( ' , 1 2 1 ' , ' , 1 2 , ' ) ~ ' » F7 . 3 ) 

END  IF 
END 


v  %  *  *  %  9f  #  %  %  %  %  %  t  %  *  %  %  %  %  %  %  %  %  %  %  %  %  %  %  *  *  *  1 1 *  *  *  #  >!<  *  *  *  *  *  *  * *  -!<  *  *  *  *  *  *  *  *  *  *  *  *****  *  *  *  * 

CUBE  OUT  I N  E  D  E 1 '  I N  E  -AIT  E  R  N  AT  I V  E  S 

THE  PURPOSE  OF  THIS  SUBROUTINE  IS  TO  DEFINE  THE  VARIOUS 
ALTERNATIVE  LOCATIONS  FOR  ENDPOINTS  J.  AND  2.  THE  PROGRAM 
IS  PASSED  THE  COORD  INATES  OF  THE  ENDPOINTS  AND  THE  STFP-S1  7M 
AMD  OUTPUTS  A  TWO  BY  TWEN  T  Y -SEVEN  ARRAY  CONTAINING  I  HE 


I V  I:  .  I  II 

n  i ■  si/i:  i:n  ins i ani  i 


1.  (  ll!|;  ;!l  ;  i  . |  j.  ■;  |l!  i  Mi  r',1  i  !  I  If !A  ;  :  ',,'L 

l:  on  a  ! ; i  i i : ;  mi  is  :'ii  cm  mm  ami  ;  \  s.i  i  1,1  r 

r  WITH  Till  ENDED  I  N  I  A  I  I  IS  Cl  NILE. 

(.; 

C 

C  #  #  %  %  # :{'  >!'.  %  % H :  r  I:  ;l  >1 :  %  X 1 # 1 *  X  *  #  *  #  #  t  *  *  #  X  X  X-  He  He  t  X  He  He  H  He  He  1 1 #  He  He  He  He  He  t  He  He  He  He  He  X-  He  H<  He  He  He  He  X- 1  He  *  He  >!■ 
C 

C  VAR  I  AM  i:  Dl  I  INTT  I  UNB 

C 


r 

REAI 

A 1 :  II  ABI  E  S 

♦ 

c 

FNDX < 1 

) 

HIE  X 

COORDINATE  Of 

ENDPOINT 

I 

c 

ENUY  (  .1 

) 

l  Hi  Y 

-COORDINATE  Of 

ENDPO I NT 

:i 

c 

END 3 ( 1 

) 

HIE  Z 

-COORDINATE  OF 

ENDED I NT 

T 

c 

SIZE- 

THE 

STEP 

-S IZE 

c: 

AL.TX  (  :i 

v  J ) 

^  THE 

X-  COORD  I  NATE 

OF  THE 

jri-i 

AL  TERN  AT  1VL 

FOR 

ENDED I NT 

c 

Al  l  Y  (  I 

j,  J  ) 

“  THE 

Y -COORD  I  NATE 

OF  THE 

jin 

ALTERNAT I  OF 

FOR 

ENDPOINT 

c 

ALTZ  <  I 

y.l) 

"  THE 

Z ■ COORDINATE 

OF  THE 

JII-I 

ALTERNATIVE 

FOR 

ENDPOINT 

C 

c 
c 
r; 
c 
c 
c 
c 
c 
c 

C*##*####**##*##*####***#**#*#*#***###*##«##**###*####*###*#*#***##**# 

C 


INTEGER  OARInM.ES: 

ENDPT J. -  THE  NUMBER  CORRESPONDING  TO  ENDPOINT  1. 

ENDPT 2=  THE  NUMBER  CORRESPONDING  TO  ENDPOINT  2 
P ■■=  THE  INDEX  FOR  THE  ALTERNATE  LOCATIONS 

C  H  A  R  A  C  T  E  R  V  A  R  T  A  B  l„  E  S  : 

SI  AIR T ( I ) “  THE  STATUS  OF  POINT  If  EITHER  FIXED  OR  UNFIXED 


S  U  B  R  OUT  I N  E  D  E  F  A  I...  T  <  E  N  D  X  E  N  D  Y  E  N  D  Z ,  S 1 7.  E  ,  E  M  D  P  T I  ,  E  N  D  p T2  ,  S  T  A  T  P  T  ,  A  L  T  X  ,■ 

S  AI..TY  v  ALTZ  ) 

R E  A  I.  E  N D X  <  2 )  r  E N D  Y  <  2  )  ,  E  M D Z  (  2 )  ,  S I Z E  ,  A  L  T  X  <  2  ,  2  7 )  ,  AETY  <  2 , 2 7 ) 

X  ALTZ ( 2  , 27 ) 

I N  T  E  0  E  R  E  N  D  P  T 1  ,  E  N  D  P  T  2  ,  P 
CII A  R  A  C  T  E  R  S  T  A  T  P  T  (SO)  *  I 

C  THESE  LOOPS  LOAD  THE  ALTERNATE  LOCATIONS  FOR  ENDPOINTS  1  AND  2 

DO  SCO  T::'  .L  1 2 
P-0 

DU  310  S=-l,l 

DO  320  R"-.l ,  1. 

DO  330  Q  -  1 ,  .1 
P  “PTL 

A I  r  X  C  T  v  1 1 )  E  N  D  X  (  I  )  ~  S  *  S I Z  E 
A  I...  T  Y  *  T  ,  P )  ~  E  N  D  Y  (  T  )  ••••  R  #  S 1 2  E 
A  L  I Z  (  T  ,  I-' )  -  L.  N  D  Z  (  T  )  -  0  *  S I Z  E 
330  CONTINUE 

320  CONTINUE 

310  CO NT f NUE 

300  CONTINUE 

C  If  ONE  01  THE  ENDPOINTS  IS  FIXED,  THEN  IT  HAS  NO  ALTERNATE  POSITIONS 
IF  <  S  I  A  1 1  •  I  <  l  NDI '  I  1 )  » 1 : 0  .  '  F  '  >  THEN 


Hi  I  T  I  ,  H  )  I:  Nil  T  \  I 

A!..  17.  (  !  ,  G  ):“(•■: Nil  Z  •;  I  ) 

COM!  I  rilll. 

I- NH  U 

:i  i  (  st  AT!  *i  •:  i  hUn  'i >  j-;o,  ■  i  •'  >  hi  in 

un  a ;  i  o  s  v  ?.  / 

ni  i  m  ?,'•) -  :  mux < 2  ■ 

A!,  i  Y  I , : > 1  i  NUV  (  2  > 

ALT/.  (  2  ,  )  -  I  NTH  (  2  > 

CONI  TNi.ll. 
i- nd  i  r 
END 


C  X  X  X  X  X  X  X  X  X  X  X  X  >!<  *  X  X X  *  >k  :*<  X  X  X  %  %  %  %  it:  %  %  *  *  X  X  X  :{<  :{<  x  X  X  X  t  X  X  X  X  X  X  X  X  X  t  X  X  X  X  X  %  *  *  X  %  X  X  X  X X  X  % 

C 

C  0  LI  B  R  0  i  .l  T I M  I-  01 1 1. 0  K  -  A I .  I  E  R  N  A  T I  y  E  B 

C 

C  THIS  SUBROUTINE  CHECKS  THE  VARIOUS  ALTERNATIVES  TO  SEE 

C  IF  THEY  LIE  IN  THE  REGION,  THEN  TESTS  THE  REMAINING  GAMUT  OATES 

C  TO  SEE  II-  THEY  INCREASE  OR  DECREASE  THE  DISTANCE  BETWEEN  THE 

C  ENDPOINTS.  IF  THE  DISTANCE  IS  INCREASED,  THIS  PROGRAM  CALLS 

C  SUBROUTINE  ' FIND  -DISTANCE "  TO  SEE  IF  THE  OVERALL  MINIMUM 

C  DISTANCE  IS  INCREASED  OR  DECREASED.  IF  SO,  THE  ENDPOINTS 

C  ARE  MOVED  TO  THE  ALTERNATE  LOCATION  THAT  MOST  INCREASES  THE 

C  OVERALL  MINIMUM  DISTANCE  WITH  TIES  GOING  TO  THE  ALTERNATIVE 


THAT  INCREASES  THE 
THOSE  NEW  LOCATIONS 


ENDPOINT  DISTANCE  THE  MOST 
TO  THE  MAIN  PROGRAM, 


IT  THEN  PASSES 


%  XXX  X-  t  X  X  xc##*#*######  X  X X X  :{<  *  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  x  xxxxxxxxxxxxxxxxxxxxxxxxxxx 
:  VARIABLE  DEFINITIONS 


REAL  VAR I i 
X(I>“  THE 
Y  <  I )  THE 
ZCt)“  THE 
CARD-  T Hi: 
AI..TX  ( I ,  J )  " 
ALT  Y  (  I , J )  ■ 
AL.TZ  ( I  ,  J )  • 


X  -CO 0 R 0 IN A T E  OF  POINT  I 
Y -COORD I NATE  OF  POINT  I 
2.  ■COORDINATE  OF  POINT  I 
M I N I M IJ  M  D I S  T  A  N  C  E  B  E  I  W  E  E  N  A I .  I. 
:  THE  X -COORD  I  MATE  OF  THE  J  l  l 
THE  Y  COORD  I  NAIL  OF  Tl  IF  J  TI¬ 
THE  7  -COORD  I  MATE  OF  (HE  Jll 


U..L  PAIRS  OF  (UNFIXED)  POINT! 
ITII  ALTERNATIVE  FOR  ENDPOINT 
I T  H  A I .  T  E R N  A  T I V E  F 0  R  E N  D  P 0 1 NT 


1 1  IE  .Mil  ALTERNATIVE 


ENDPOINT  I 


NEWD~  Ti  ll  NEW  MINIMUM  DISTANCE  BE  I  WEEN  ALL. 


THE  POINTS 


Cl IKX ( I ) “  THE  X -CO ORDINATE  OF  POINT  I  BEING  CHECKED  AS  A  MEW 
LOCATION  SCHEME  BY  SUBRUUI  INK  '  FIND- DISTANCE  ■' 

C  FI  K  Y  (  I )  =•=  1 1 1 E  Y  -  C  0  0  R  D  \  N  A  T  E  0  F  P  0 1 N  I  I  B  E  \  N  G  C  FI  E  C  K  E  D 
CHKZ(I)-  THE  Z- COORD I NAT F  OF  POINT  I  BEING  CHECKED 
DUM  J. ,  DU  M2 ,  DUM.I~  DUMMY  VARIABLES  THAT  ARE  RETURNED  BY  SUBROUTINE 

MIND  -DISTANCE"  BUT  ARE  NOT  USED  IN  THIS  SUBROUT  INI! 
MIND  THE  MINIMUM  DISTANCE  CALCULATED  BY  SUBROUTINE  " FIND  -DISTANCE ' 


M  I  ND 


I  (H-:  ill!:  i  .  0  C  A  I  I  (It/  '..CII!  i'ii  nl. !  lNi.il  1  I.IIK).  i-Clll.i  ,  ..i'H 

MAXIMD-  I  III  VAR  I  Aid  E  USED  m3  A  I  1 1.  BREAM  R  I  01:  AL  I  I  RNATIVE.S  '-J  ill! 

i.  n  e  n'i  i  c  a  i  m  i  n  :i  h  u  h  d  i  s  r  a  n  c  r  s 

E  Nil  PI  IP-  THE  DISTANCE  BETWEEN  ENJlPIJINTS  J.  AND  2 

YE; ED  (  [  ,.  J  >  ■  THE  VALUE  01  A!  IX  <  I  ,• . I  >  I  N  RED  GRIT  N -BLUE  SPACE 

Y  G  R !:.  E  N  (  I  ,  J  >  Til  i:  V  A I ..  (J  E  0 1  A  I...  T  Y  (  I ,  J  )  I N  R  0  D  S  P  A  C  E 

Y  D I  U  I  -  ( I  f  J  )  -  III  E  V  A I .  U I  0 1-  A  L  I  2  (  I  )  I N  R  G  B  3  P  A  C  If: 

INTEGER  VARIABLES? 

E N D P T  I  ~  1 1 1 E  M U M HE R  A 3 1 ! 0 C I A  I  E D  W  f  I  El  E N D P 0  I N  T  i 
ENDI  'T2:::  THE  N  UNDER  ASSOCIATED  WITH  ENDPOINT  2 
N THE  N UNDER  OF  I  T) I  NTS 

DUM4  y  DUNS""  NOT  USED  IN  T  HIS  SUBROUTINE 
C  H  A  R  A  C  T  E  R  V  A  R I A  B  L  E  S  : 

ST A TPT ( I > “  THE  ST  ATUS  OF  POINT  I ?  EITHER  FIXED  OR  UN!  IX I  D 
LOCALT  ( I  ,  J )  •"=  THE  LOCATION  OF  THE  JTH  ALTERNATIVE  FOR  ENDPOINT  I  5 
EITHER  IN  OR  OUT 

DUMA-  A  VARIABLE  REQUIRED  FOR  SUBROUTINE  '  RGBOUT '  BUT  NOT  USED 
IN  THIS  SUBROUTINE 


t  *  *  :{<  *  *  *  #  *  #  %  %  %  *  t  *  *  *  :L  *  #  #  #  *  *  %  %  *  *  t  *  #  #  %  %  %  %  #  #  *  X  %  %%%%%%%% 

C 
C 

SUBROUT  I NE  CHKAI..T  ( X  ,  Y  ,  Z  ,  N ,  CARD ,  AI..TX  ,  AL  T  Y  ,  ALTZ  ,  STATPT ,  ENDPT  .1. 

X  ENDPT2, NEWD ) 

REAL  X  <  50 ) , Y ( 50  > , Z  <  50  > ,  CAPD ALTX  ( 2  27 )  ,  ALTY  (  2  27 )  y  ALTZ  ( 2  ,  27 )  y 
X  NEWD  y  CHKX  ( 50  )  ,  CHKY  (  50 )  y  CHKZ  (  50 )  y  BUM  1  (  2 )  y  0IJM2  <  2 )  »  BUMS  (  2  )  r 

X  MI N D  y  M A X E N D  ,  E N B P T D  r  Y R E D  (  2 , 27)  » YG R E E N  (  2  ,  2 7  )  ,  Y B L U E  (  2  ,  2 7  ) 

I  NT  EGER  ENDPT  1  y  ENDPT  2  y  N  ?  DIJM  4  ,  BUMS 

CHARACTER  STATPT  <  50 )  >Y- 1 ,  L  DC  ALT  C  2 27 )  #3  y  Ji  U  M  6 1 A 

NE WD-CAPD 

MAXEND-CAPD 

DO  400  I ~ 1 y M 

CITKX  <  I )  =  X  <  I ) 

CHKY <  T )=Y( I ) 

CHKZ ( I ) -Z ( I ) 

400  CONTINUE 

T  THIS  LOOP  CHECKS  ALL  54  ALTERNATIVES  TO  SEE  IE  THEY  ARE  IN  OR  OUT 
:  OF  THE  REGION 

DO  4  .10  S- 1 ,2 

DO  420  T-.1.  ,27 

CALL  RGBOUT (A!.TX( S ,•  T  )  r  ALTY  (  S  T  )  ,  ALTZ  < S  ,  I  )  1  ,  YRI.  D  <  S  T  )  , 
X  YGREEN ( S , T ) , YBLUE ( S , T ) , DUMA ) 

C  A  L  I...  C  H  E  K  P  T  (  Y  R  E  B  (S,T),Y  G  R  E  E  N  <  S  y  T  )  ,  Y  B I .  U  E  <  S  ,  I  ) 

X  LOCALT ( S , T ) ) 

120  CONTINUE 

HO  CONTINUE 

DO  430  T - I ,27 

IF  (LOCALT (1,T)  » E  0 .  'IN')  THEN 
DO  440  S=l,27 


.1.1  1‘ U  \  ;  !  !  r-'.'i  ( 

HI:  1  W!  I  N  TULA 


n  THAT  t;i;i  ;;tam(.t 


K  I  .i.  I  1  . 


Ill'  N  i  nl  i  III  nil 


n  i .  u  bj:h  i  anc 


I  fin:  1  l.i  -  URT  (  (  AL  I  X  C  1  ,  1  >  hi  TX  (  2  ,  S  )  )  **2} 
(  ALT  Y  (  1  ,  I  )  AI.TY  <  2  L. )  )  *  * 2  }•  <  ALT?  (  :l  ,  T  )  ■■■ 


(•HI  ATI  I:,.  Till  i!  CHI  CK  TUI.  OVLRAI.I 

rr  (cm  up  i  n  hit,  ml  kin  >  tit  cm 


IN J MUM  JLI 


Cl  IK 

X  (  E 

i 'DHt  :i  > 

-  AL 

1  X  (  1  V 

T  ) 

Cl  IK 

y  •:  l 

nuut  .1.  ) 

- :  A  l 

T  Y  (  1.  y 

T) 

Cl  IK 

2  ( E 

mu'll  > 

-  AL 

T2.(  1 

T  ) 

CHKX ( F 

NDPl 2 ) 

"AL 

IX  ( 2 , 

S ) 

Cl- IK 

Y  ( E 

NBPT2) 

::::  fi  | 

TY  (  2  y 

S ) 

Cl  IK 

Z  <  E 

NDP f  2  > 

"At. 

!  7.  (  2  y 

S) 

CAL 

1.  F 

INDD ( CHKX 

;  r  C!  IKY 

r  Cl 

s  n  LJ  M 1.  v  IT  u  H  2  ,  D  LJ  M  3 ,  D  U  M  4  ,  BUMS  r  IT  I  !  M  6  ) 

I  I-  THE  MINIMUM  DISTANCE  I S  Tin:  SAME  but  with  a  greater  endpoint 
DISTANCE  P  OR  IF  THE  MINIMUM  DISTANCE  IS  BETTER  f  THEN  HOOF  T MI¬ 
EN  DP  0  I  NTS 

IF  ((MIND  .GT .  NEWD)  .OR.  (MIND  .EQ. 

X  NEWD  .AND.  ENDPTD  .GT.  MAXEND))  THEN 

NEWD --MIND 

MAXEND ^ENDPTD 

X ( ENDPT 1 ) ~ALTX ( 1 » T ) 

Y ( ENDPT 1 ) "ALTY ( 1  r  T  > 

Z(  ENDPT  l.)~ALTZ(l,T) 

X  ( E.NDPT2 )  "AI..TX  ( 2  y  S ) 

Y ( EM DP  12 ) "AL IY ( 2 , S ) 

Z ( ENDPT2 ) ~ALTZ ( 2  y  S  > 

END  IF 


END  IF 


END  IF 
CONTINUE 


END  IF 
CONTINUE 
END 


C  *****  *  *  * 1 1 1 *  *  ***  *  t  *  #  *  *  t  *  %  *  %  *  %  %  %  :{<  *  %  % :«  %  t  * *  *  ***  *  *  *  %  %  %  %  ******  *  *  *  *  *  *  *  *  *  *  *  * 
D 

c  S IJ  B  ft  0  U  r  I N  E  C  H  E  C  K  -  R  0  I.  N  I 

:  THIS  SUBROUTINE  CHECKS  POINTS  TO  SEE  IF  THEY  ARE  IN  OR 

OUT  OF  THE  REGION.  IT  CHECKS  THEM  AGAINST  A  SET  OF  EQUATIONS 
"  (CONSTRAINTS)  OF  THE  FORM  AX  T  BY  +  CZ  <  K  WHICH  MUST  BE 

:  SUPPLIED  BY  THE  USER f  AND  THEM  OUTPUTS  A  CHARACTER  VARIABLE 

■  WHICH  IDENTIFIES  WHETHER  THE  POINT  IS  'IN'  OR  'OUT.' 


note:  HUS  subroutine  is  for  linear  constraints  only,  non 

I  INEAR  CONSTRAINTS  MAY  BE  USED  BY  CHANGING  THE  EQUATION  FOR 
V  <  S ) v  BUT  THEY  MUST  BE  CONVEX  IN  ORDER  FOR  THE  ALGORITHM  TO 


(/  *  %  if.  t  #  >}: »  *  *  *  *  *  )j:  #  >!■:  j}:  *  #  *  if.  if.  #*:#  *  #  *  #  *  #  *  #  *  #  *  *  %  if  #***  *  *  *  *  *  *  *  *  *  *  *#*##*##  *  ##*###  *  * 
(' 

C  VAR 1 A DU  DL  L  IN  I  T  KINS 


r 

PI  A!  V 

aim  ALi  i-:s: 

( 

XI  UORB 

III!  X  C 

OORB I NA  T  L 

01  1  III-  POINT  BEING  TL  ST  ED 

y  pool :  i  ’ 

1 1  it;  i  i 

OORBINATL 

OF  THE  POINT  BEING  TESTED 

C 

7 CO OR H 

I  HI  1 

OORB I NATL 

01  THE  POINT  BEING  TESTED 

i: 

A  <  I  )  • 

Tin  nil  i  r 

ICE  ENT  OF 

X  IN  CONSTRAINT  I 

i: 

It  <  |  ) 

I  IIP  (ONI 

I  CLIENT  OF 

Y  IN  CONSTRAINT  I 

c 

i :  (  i  --  •• 

hi:  i:olfe; 

1.0  IFN  I  or 

7.  IN  CONSTRAINT  I 

C 

K  <  I  )  - 

THE  V  A 1  01 

01  THE  R 

I GILT  HAND  SIDE  IN  CONSTRAINT  I 

c: 

v  (  r ) 

THE  VALUE 

01  CONSTRAINT  l  EVALUATED  AT  THE  POINT 

c 

BEING  TLS 

1  I.  IT 

L 

c 

INTEGE 

R  VAR  I ABI 

ES? 

I-™  THE  HUM  HER  HE  FACES  (CONSTRAINTS)  ON  THE  REGION 
CHARACTER  VARIABLES ? 

POINT"-  THE  LOCATION  OF  THE  POINT?  EITHER  IN  OR  OUT 


C 

C 

SUBROUT INE  CHEKPT  <  XCOORB YCOORB  >•  ZCOORD ,  POINT ) 

R  E  A  L.  X  i :  0  0 1  •{  B ,  Y  C  0  0  R  D ,  Z  C  0  0  R  B ,  A  ( 1 0 ) ,  B  ( 1 0 )  ,  C  ( 1 0 )  ,  K  <1 0  )  ,  V  (  1 0 ) 

INTEGER  F 
CHAR AC  I ER  POINTER 
C 

C  •—  IMPORTANT  NOTE  TO  USER 

C 

C  THE  USER  MUST  SUPPLY  THE  NUMBER  OF  FACES  FOR  EACH  NEW  PROBLEM  HERE 


C  THIS  I. OOP  INITIALIZES  THE  COEFFICIENTS  AT  ZERO 
BO  500  0 " 1  ,  F 
A  (  S )  •"  0 . 0 
B  (  S  )  ~  O  .  0 

r:  ( s )  -o o 

K  (  S  )•().() 

V  (  s  )  "0  0 
500  CONTINUE. 


•  -  IMPORTANT  NOTE  TO  USER  — 

T HI  USER  MUST  SUPPLY  THE  LOCATIONS  FOR  THE  CONSTRAINTS  IN  THE 
TORN  AX  I  BY  f  CZ  P  K  I . E .  HE  MUST  DETERMINE  AND  PROVIDE 
THE  VALUES  01  A(S),  IKS),  CIS),  AND  K(S)  FOR  S"  1  TO  F .  FOR 


ltj  in 


j ,<s  i  r;'i  v  .  i  niK  :  lpst  c  i 's  i 1  :A !  n  i  is  ox  i  :> ;  i  4.  _  s.  ii.in 

f i*  i  > .  • '  .  I M'  |  t :  (  .1  » :  /;  ,.  .11  K  (  :l  I!  'll;!  SI  OGL'D  ('(INST  Hr1.  IN  I 

TS  X  X  1  .  THEN  A  ( 2. )  I  y  B  (  ?  >  0  v  (",(:•»  )™0,  AND  K  (  2 >  .!  .  I  F  A 
CONST  RAJ  NT  HAS  A  ORi  AIER  THAN  SIGN,.  Till:  N  REVERSE  Till:  SIGN  BY 
N II G A  T  l  M f  ,  T 1 1 F.  C 0 F F  F  J  I  I F t) T S  ,  !  0 R  1  N S  I  A M ( T !.  ,  X  I  Y  T  7  2  7  B F SO M F S 

•  X  ■  Y  -v  <  3.  SO  A  •'  3 )  ••  •  •  I  BIS'  I  >•  C  <  3  ) :  - -  1  >•  ANU  K  <  3  )  1,  II- 

THF  CONSTRAINT  IS  AN  EQUAL  I  I'Y  CONSTRAINT.  I HFN  FORM  I  HO  MOD  AT  I  ON 
WITH  ]  NFC'LIALIT  ICS.  I  OR  INSTANCE.  Y  I  7  -  G  BECOME  S  Y  !  Z  _  S 

A  N  l.t  -  Y  7  2  ■  S  ,  S  0  A  (  4  >  =  0 ,  B  (  A  >  -  - 1  .  C  (  4  )  -  I  ,  K  (  4  )  =  5 .  A  ( 5 )  0  >■ 

B<S>  '  I.  0  <  F, ) AND  K  ( ■  5.  IF  YOU  KAMI.  ANY  PROBLEMS, 

C  ON  SO  FT  YOUR  LOCAL.  MATH  EXPERT. 


A  ( 1  ) 

=  1 

* ) 

B  <  2 ) 

- 1  . 

0 

C  <  3 ) 

I  . 

0 

A  (4) 

”  —  .1. 

.0 

B  (  5 ) 

=  ■  .1 

.  0 

C  <  6 ) 

=  -•.1 

.  0 

A  ( 7  > 

•~™1 

.0 

B  ( 7 ) 

.0 

C  <  7 ) 

.0 

K(  :t ) 

“50 

.0 

K  ( 2 ) 

=  160.0 

K  <  3 ) 

=20 

.0 

K  <  4 ) 

=0. 

0 

K  <  5 ) 

=0  . 

0 

K<6> 

=  0. 

0 

K  ( 7 ) 

=5-2 

.3 

C  THE  CONSTRAINTS  AS  THEY  ARE  SHOWN  ABOVE  CORRESPOND  TO  CARTER  AND 
C  CARTER'S  CONCEPT  OF  THE  FEASIBLE  REGION  FOR  THE  COLOR  SPACING 
C  PROBLEM  AND  LOOK  LIKE  THIS  IN  RE By  GREEN,  ANU  BLUE  COORDINATES: 

C  RED  X  SO 

C  GREEN  <  .1.60 

C  BLUE  <L  20 

C  RED  >  0 

C  GREEN  >  0 

C  BLUE  >  0 

C  RED  }■  GREEN  E  BLUE  >  2.3 

POINT-' IN' 

C  THIS  LOOP  CHECKS  THE  POINT  AGAINST  EACH  OF  THE  CONSTRAINTS 
BO  510  S~1 , F 

V  <  3 )  --  A  ( S )  5CXC00RD-EB  <  S )  tfYCOORD  1C  ( S  )  rZCOORD 
IF  <v-:s>  .GT,  K  ( S )  )  THEN 
POINT  '01.1  I  ' 

GO  TO  520 

ENT  i  || 

1.0  CONTINUE 
20  END 


C  *  %  *  t %  %  %  %  %  *  t  %  *  %  *  *  y  *  *  #  *  *  :fc  #  %  >K  $  >K  *  *.  t  H<  t  $  *  *  Y 1 * Y  Y  %  t  >1 *  $  X>:  *  *  *  >1  #  *  T<  *  *  *  *  *  *  >(< !)' 


cudpuu  i  i ! 1  j 


' . ( l * M  : ( ! 1 1  !  !  !  '|  !  :  ;•  !  I  :  ! 

mi  i M.1K1  ! i': ,i-  m  in i s  s.u  brgiu  i. {'!(•.  i s  to  i  iv  points  and 
m:  ii  kmjmi  u ;  i  r :  t  1 1 1'  i :  or  hoi  thi;  program  is  hone,  the  logic 

ASSOC  I  AT  I.Ti  0.1(1)  I  I  I  Of,  l-M  IN  |  G  |  0  THAI,  ONCE  THE  MINIMUM 
DISTANCE  IS  MAXIMIZED,  THAT  HI  STANCE  JS  MO  LONGER  UNDER  CONSIDER 
AT  ION ,  SO  ONE  01  1 1 II  !  NDP0IN1S  IS  I  I XI  D  AT  ITS  LOCATION  AND  THI 

SECOND  CLOSES  I  PAIR  01  'TUNIS  IS  MAXIMIZED.  IMIS  CONTINUES 
UNI  I!  All  POINTS  ARE  11X10,  AT  WHICH  TIME  THEY  ARE  UNFIXED 
TO  SEE  !  i  ANY  'IMPROVEMENTS  CAN  DC  MADE.  THE  PROGRAM  STOPS 
WHEN  NONE  CAN  BE  MADE. 


C  *  *  *  *  *  :  *  *  *  IT  *  *  *  * * *  *  *  *  *  *  *  *  *  *  *  #  *  *  *  *  #  :{<  * :)( *  *  *  ***  *  *  ******  *  *  *  *  *  *  '!•:  *  *  *  *  :{<:  *  *  *  *  *  *  *  * 
C 

C  0  A  R  T  A  D  I...  E  D  E I  *  I N  I  T  1 0  N  S 

C 

c  real  oariadi.es: 

C  X  (  '[)-•  THE  X -COORDINATE  OF  POINT  T. 

C  Y  <  J  >  "■  THE  Y -COORD  I  NATE  OF  POINT  I 

C  Z(I>=  THE  Z -COORD I NATE  OF  POINT  I 

C  D  U  M  .1.  I)  I.J  M  2  ,  D  U  M  3 ,  D IJ  M  A  ~  HUM  M  Y  0  A  R I A  B  I...  E:  S  R  E  Q IJ  I.  R  E II  F  0  R  SI  J  B  R  0 1 J  T  I.  N  E  '  F I 

C  D I  STANCE"  BUT  NOT  USED  IN  THIS  SUBROUTINE 

C  ME  MX  <  I. )  “  THE  X- COORD  I  NATE  OF  POINT  I  FROM  THE  PREVIOUS  ITERATION 

C  MEMY < I ) “  THE  Y -COORDINATE  OF  POINT  I  FROM  THE  PREVIOUS  ITERATION 

C  ME HZ < I ) =  THE  Z- COORD I NATE  OF  POINT  I  FROM  THE  PREVIOUS  ITERATION 

C  SIZE-  THE  STEP -SIZE 

C 

C  INTEGER  VARIABLES: 

C  N~  THE  NUMBER  OF  POINTS 

C  ENDPT.1  ~  THE  NUMBER  ASSOCIATED  WITH  ENDPOINT  I 

C  ENDPT2-  THE  NUMBER  ASSOCIATED  WITH  ENDPOINT  2 

C  NUMFIX-  THE  NUMBER  OF  POINTS  THAT  HAVE  BEEN  FIXED 

C 

C  C  l-l  A  R  ACT  F.  R  V  A  R I A  B  I..  E  S  : 

C  GTATPT ( I ) ~  THE  STATUS  OF  POINT  I?  EITHER  FIXED  OR  UNFIX EH 

C  STATUS-  THE  STATUS  OF  THE  PROGRAM;  EITHER  DONE  OR  RUM 

C 
C 

C  *  *  *  *  *  *  *  *  *  *  ******  *  *  *  *  *  *  *  *  *  *  *  *  *  ******  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *  *•  *  * 

c 

c 

SUB R OUT  I N E  F I X  P T  <  X  ,  Y  v  7. ,  N ,  S T A T F ' T  ,  E N D i  'T  :l  ,  E N 0 P T 2 ,  S  I'  Z E  ,  S T  A T U S  ,  H T  M  X  , 

8  MEMY, ME HZ) 

REAL  X ( 50 ) , Y  < 50 ) ,Z ( 50 ) , DUN  I  , BUM 2 ( 2  > . BUM 3 < 2 ) , BUM 4 ( 2 ) , ME MX ( 50 ) , 

8  M  E  M  Y  (  5  0 )  ,  M  E  M  Z  <  5  0  )  ,  S I  /  L 

I N  T  P  0  E  R  N  ,  E  N  D  P  T 1 ,  E  N  D  P  T  2  ,  N  IJ  Ml"  I X 
G  H  AR  A  C  T  E  R  S  T  A  T  P  T  (  5  0 )  *  1 ,  S  T  A  I  IJ  S  *  A 

C  —  IMPORTANT  NOTE  TO  USER  — - 


IF  THE  USER  WISHES  TO  CHANGE  THE  MAXIMUM  STEP-SIZE. 


E  MUSI  MAKE 


I  TMi  p.H.V-.'i:!  ill  Rj  AS  Ui-.l  I  AS  J  N  Kli  LPiiN  I 'EUOROi  i . 

(., 

S 1 7.  E  =••  J.  'i  0 . 0 
NU HI  IX-  0 

c  mis  i. oop  coon  is  ri-i::  number  or  points  that  are  al  ready  fixed 

DO  600  J  1  ,  i? 

T !  S  C  A  T  PT  (  I  )  ,  I  0  ,  •'  P  •'  >  THEN 

n  1 1 m :  -  j:  x  -  n  1 1  it  r  i  x  :•  i 

P.ND  IF 

600  CONTINUE 

c: api  .  r  i ndd  •;  x ,  y  ,  x n  ,  btatpt  uumi  ,  bum:?  ,•  dijm3  ,•  duht  ,  enbpt :i.  ,  pnbpi  :> 

X  STATUS) 

C  THESE  NEXT  TWO  IP  STATEMENTS  FIX  A  NEW  POINT 
I F  (  ST  A  l  PT  (  EM  DPT  .1.  )  <•  EO  *  '  U  ■'  >  THEN 

NUMFIX-NUMF IX  I  1 
ST  A  LPT  (  ENDPT  X  )  “  '  F  ' 

GO  TO  610 
END  ir- 

I F  <  S  T  A  T  P  T  (  E  N  D  P  T  2  >  .  E  Q .  '  U '  >  T II E  N 
N  U  M  F I X  M  U  M  F I X  M. 

ST ATP T  ( ENDPT2 )  “ '  F  ' 

END  IF 

C  IF  ALL  THE  POINTS  ARE  FIXED,  THEN  CHECK  TO  SEE  IF  THEY  ARE  IN  THE 
C  SAME  LOCATION  AS  BEFORE 
610  IF  ( NUN FIX  .EG,  N)  THEN 
DO  620  1=1 , N 

IF  ( ME  MX  ( I )  ,  NE  ♦  X(I)  .OR.  MEMY  ( I )  .NE»  Yd)  .OR. 

X  MEMZ < I )  .NE.  Z ( I ) )  THEN 

STATUS™ ' RUN ' 

00  TO  630 

ELSE 

STATUS  - ' DONE ' 

END  IF 

620  CONTINUE 

C  IF  THE  POINTS  HAVE  MOOED  SINCE  THE  LAST  ITERATION  ,  THEN  UNFIX  A I  I. 
C  THE  POINTS,  STORE  THEIR  LOCATIONS  IN  MEMORY,  AND  RIJN  THROUGH  THE 
C  PROGRAM  AGAIN 
630  DO  6 TO  I--.1,N 

STATPT  ( I  )  '  U  •' 

ME MX ( I ) ~X ( I ) 

MEMY ( I ) =Y < I ) 

MEMZ < I ) =Z ( I ) 

6  TO  CONTINUE 

END  IE 
EN  0 


C  :K  1 1***  %  t Z  *  >t-M  %  %  %  %  *  :f<  *  %%  %  X  ##*)((  %  %%  X  ZXZZ  #  **  Y  * 

C 

C  SUBROUTINE  PRINT 

C 

C  THIS  SUBROUTINE  PRINTS  Till;  LOCATION  01  AIL  N  POINTS  TO  1  HI 


it  ■  •  i ■  i  mi i ;  '  i 


om  '  ■  . 


i  i  i  i  •  ■  , . 1  ■  :  (i  si  •!.:(  ii  i  T  '  ■:  .1  i  !  K!  :  1 1  i  ni  I  1  \  T  i 

r, 

c 

(.:*  *  X  %  %  X  X  X- X X  t  X  •;  *  X  X  X  t  X  X  X  X  t  X  X  X  tX  X  X  X  t  X  %  t  X  X  X  X  X  X  X  X  X  X  X  X  X  X  %  X  txx  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X 

i; 

C  V  ft  I'M  A  HI  I:!  Til  I  I  N  [  I  I  (IK'S 


.:  ri.ai 

I  A! :  1  AT: 

1.  I'M 

3  X  1  :• 

Tin  X 

-COORD IN A  I  E 

OF 

PO  I  NT 

:  y  <  ;i  >  ^ 

mi  y 

- GOURD  1  NATE 

or 

pc  i  :i  nt 

z  ( r  > 

i  hi:  7. 

COORDINATE 

01 

POINT 

:  IN  IT  G 

I:  VAR 

l A DEES  l 

c 


N  rill;,  i'll  111  ITS  R  OF"  POINTS 


c  X  X  X  X  X  X  X-  X  X  X  X  X  X  X  X  He  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X  X-  X  X  X  X  X  X-  X-  X :  *  X  X  X  *  x  X  X-  X  X  XXXX  X  X  X  X  X  X  X :  X  X  X  X  X 

c 


S  L  J  IT  FT  0 IJ  T I N  E  P  R I  N  T  <  X  ,  Y  ,  7. ,  N ) 

REAL  X  ( 50 > , Y  <  50 )  ,  7.  <  50 ) 

INTEGER  N 
WRITE <3, 700) 

DO  71.0  T  =  :l.  ,N 

WRITE <3,720) I ,X< I ) ,Y< I )  ,Z< I > 

CONTINUE 
WRITE <3 ,730) 

FORMAT  (  '  POINT  '  ,  ?X  ,  '  L# '  ,  13X ,  7  IJ*  '  ,  13X  ,  '  V*  '  ) 
FORMAT ( 3X , 1 2 , 7X , F8 . 3 , 7X , F3 . 3 , 7X , F8 . 3 ) 

FORMAT ('  ' ) 

END 


710 

700 
720 
730 

C 

c 

( J  X  X  X  X  X  X  X  X  X  X  *  X  X  X-  X  XX  X  **  X  X  X  X  XX  X  *  X  >}c  X<  #  *  X  X  X  XX  X-  X  X  *  X  XX  X- X  >l<  *  *  *  X  X  *  X  X  X  *  X  X  ft#  X-  *  X  X  X  *  * 

c 

C  S  U  B  R  0  IJ  T I N  E  R  (3  B  0  U  T 

C 
0 
C 
C 
C 
C 
0 
0 
C 

c 
c 
c 

cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

c 

C  VARIABLE  DEFINITIONS 

C 


THIS  SUBROUTINE  CONVERTS  POINTS  FROM  THEIR  I..*,  U*,.  AND  V* 
COORDINATES  TO  THEIR  TR I STIMULUS  VALUES  AND  THEN  TO  THEIR 
PHOSPHOR  LUMINANCES  AND  PRINTS  THESE  VALUES  TO  THE  FILE  NAMED 
' SRC OUT '  ( A . K  >  ft ♦  SPACE-OUT).  THE  USER  MAY  CHOOSE  TO  ESTABLISH 
DIFFERENT  VALUES  OF  MAXIMUM  LUMINANCE,  U -PRIME  HOT,  V  PRIME  NOV, 
AMD  CIIROMATICITY  COORDINATES  FOR  THE  GUMS  DEPENDING  ON  ITIS  OR 
HER  PARTICULAR  APPLICATION.  ALSO,.  FOR  SPACING  IN  A  SYS  I F  M 
OTHER  THAN  THE  CIE  L*U*V*  SYSTEM,  DIFFERENT  TRANSFORMATION 
EQUATIONS  MUST  BE  USED. 


68 


I. 

c 


L.S  I  Ci  ■:>,  i  I  III  v'n!  I '!  HI  I  »  I  UK  i’ll  Li  I  )  I 

1 1  S  I ;  1 1-:  (  l  >  ilil  v'n  I  .UI-.  or  |l>!:  I- OK  POINT  I 

C  VClARd  >-  HR  VAI  Ul  01  V*  FOR  POINT  I 

0  XU  HU  X  (II! '!  i  On  1  IF!  TV  COORD I  NATE  FOR  TV!!:.  KID  BUN 

l  YU  I HF  Y  CHI  Mini  l('l  1  Y  COORDINATE  FOR  T!  IF  K!  H  OIJN 

I  XT  III!  X  CHROMA  III  MY  00  OKIt  I  NAIF.  FOR  Till.  ORPIN  BON 

0  VO  ’Till.  I  ll'i.'lrinl  ill  I  Y  COORD  I  MA  IE  FOR  I  IIP  ORION  GUN 

0  XT-:-  |  III.  X  CHROMA  I!  Cl  I Y  COORD  I  NAT  I  FOR  Till  Rl  01  BON 

0  YD  I  Hl.  Y  CHROMA  I  T.C  I  i  Y  COORD  I NA 1 F  !  OR  I  HF.  Ul  II!  BUM 

C  YiH.il  Till:  MAXIMUM  !  MM  1  HANOI .  API  00  ADI  I: 

!  0 I  'M  NOT  1  I  II  .  V  A  I  .OF  OF  0  PR  [Ml  NOT 

C  l.H  'MNO  1  T  I  II:  VALUE  01  U  PRIMi:  NOT 

0  X  TOTAL.  (  i:  ) :!i?  1 1  IE  X  COORDINAIF.  OF  THE  TRISTIMUL.U3  V  ALT. I  IE  FOR  POINT  I 

C  Y  T  0  1  A  I...  (  i  )  ™  TIT  E  Y  •  0  0  0  R  Ti  t  N  A  1 1.  0 1  1 1 1 F  T I  •!  I S  T 1 H  0 1...  U  S  V  A I..  0  E  F  0  R  P  0 I N  T  I 

C  7.  TOTAL  ( I )  ••=  THE.  Z  •COORDINATE  OF  THE  TR I  STIMULUS  VALUE  FOR  POINT  1 

C  UPRIMF  ( I )~  THE  VALUE  OF  IJ  -PRIME  FOR  POINT  I 

C  VPRIME<I)”  THE  VALUE  OF  V  PRIME  FOR  POINT  I 

C  SMAL.LX  ( I  )  --  1 I  I F  X  C HR 0 M  A  J I C I T  Y  C 0 0 R D I  NAT  E  0 F  P 0 1 N T  I 

C  SMALL. Y  ( I )  ~  T 1-1 E  Y  C H R 0 M A T I C I T Y  C 0 0 R D I N A T E  0 F  P 0 1 N T  I 

C  SMALL  Z(J.)  =  THE  Z  CHROMATIC  I TY  COORDINATE  OF  POINT  I 

C  YRED(  I  >“-  THE  RED  PHOSPHOR  LUMINANCE  FOR  POINT  I 

C  Y GREEN < I >  =  THE  GREEN  PHOSPHOR  LUMINANCE  FOR  POINT  I 

C  YDLUE  (I)"-  THE  BLUE  PHOSPHOR  LUMINANCE  FOR  POINT  I 

C  K 1  ~  THE  FIRST -ROW,  FIRST  •COLUMN  VALUE  OF  THE  TRANSFORMATION  MATRIX 

C  THAT  CONVERTS  TR I STIMULUS  VALUES  TO  PHOSPHOR  LUMINANCES 

C  K2=  THE  FIRST-ROW,  SECOND -COLUMN  VALUE  OF  THE  MATRIX 

C  K3::::  THE  FIRST-ROW,  THIRD -COLUMN  VALUE  OF  THE  MATRIX 

C  K  A  -■  THE  THIRD -ROW,  FIRST  -COLUMN  VALUE  OF  THE  MATRIX 

C  KS=  THE  THIRD-ROW,  SECOND-COLUMN  VALUE  OF  THE  MATRIX 

C  K6~  THE  THIRD-ROW,  THIRD-COLUMN  VALUE  OF  THE  MATRIX 

C  DPT-  THE  DETERMINANT  OF  THE  TRANSFORMATION  MATRIX 

C 

C  INTEGER  VARIABLES: 

C  N -  THE  NUMBER  OF  POINTS 

C 

C  C  I- 1 A  R  A  C  T  E  R  VAR  I A  B  I. .  E  S : 

C  STATUS-  THE  STATUS  OF  THE  PROGRAM?  EITHER  'DONE'  OR  'RUM' 

C 


C 

c  *  *  :i<  t  %  %  %  %  t  %  #  ttz  *  X  ##  %  ))<  % *  :i<  >K ^  «  t  %  >K  %  %  -!<  *  'M  *  *  *  *  %  %  l<  :*  #  *  *  #  t  *  *  *  #  *  #  *  *  **  *  # 

C 


SUB  ROD  I  INK  RGB  01  IT  <  I  STAR  ,  OS  IAI  i ,  VST  AR  ,  N ,  YRED  ,  Y  GREEN  »  YBI..UE  ,  STATUS  ) 


REAL 


L.S  TAR  <  SO )  ,  US  TAR' 


>  ,  V  S  1  A  R  <  5  0 )  ,  X  R ,  Y  R ,  X  G  ,  Y  G  ,  X  B  ,  Y  B ,  Y  N  0  T  , 


%  Vl'MNI  1 1  ,  UI-'HNO  1  ,X  TOTAL.  (SO)  ,  YTOTAL.  (  SO )  ,?TOTAL  (SO)  ,  UPRIME  (  SO  )  , 

y,  VRR I  ME  (  SO  )  ,  SMALL.X  <  SO  )  ,  SMALL. Y  ( SO  )  ,  SMAI..LZ  (  SO  )  ,  YRED  (  SO )  , 

8  YGRL.EN  (  SO  >  ,  YBI  UE  (  SO )  ,  K I ,  K2 ,  K3  ,  KT  ,  KS  ,  K6  ,  DEL 


INTEGER  N 

Cl  lARAC'T  L  I :  ST AT US #4 


UPMNO  I  :  «  IV'.:- 3/6 
VPMNO  I  ,  <569/00 
KJ  •  XIS/YIS 
K  2  =  X  G  /  Y  G 
K3-XB/YB 
K  4  <  l-XR-YD/YR 
K5:-  <  1  •  X G  YG )  /YG 
K  6  -  (  1.  -XB-YB )  /YB 

Ii  E  1  K 1  $  K  6  ••••  K  .1.  *  K  5 I K  2  *  K  4  •  l\  2  *  K  6  +  K  3  %  K  5  -  K  3  #  K  A 
DO  800  I “ 1 f  N 

Till-: BP  ARE  Till-  TRANSFORMATION  EQUATIONS 

YTOTAL  ( I  )-=((  ( LSTAR  < I  )  +1 6  >  / 1 1 6  )  **3  )  YYNOT 
UP  R I  ME  (  I )  =  <  US  TAR  ( I )  /  ( 13#L8TAR  ( I )  )  )  TUP  M  NOT 
OPR  I  ME  ( I )  -  <  OS  FAR  (I )/( 1 3*1.  STAR  CD))  +VPMNOT 
S  M  A  L  L  X  ( I )  =  9  *  U  P  R I M  E  ( I  >  /  <  6  *  LI  P  R I M  E  ( I )  - 1 6  *  0  P  Is  I M  E  <  I )  + 1 2  ) 

SMALL Y ( I ) -4*0PRIME < I > / <  6* UP RIME ( I ) -16*VPRI ME ( I )  +  12 ) 

SMALLZ  <  I  >  ~  .1  -SMALLY  <  I  )  -SMALLX  ( I ) 

XTOTAL < I ) =SMALLX ( I ) KYTOTAL < I > / SMALLY ( I ) 

ZTOTAL  <  I )  —SMALl.Z  ( I )  # YTOTAL  <  I  )  /SMALLY  <  I ) 
YRFD<I>=<XT0TAL<I>#<K6-K5)+YT0TAL<I>*<K3*K5-K2#K6>+ 

X  ZTOTAL  <  I )  #  <  K2--K3 )  )  /DET 

YGREEN  ( I  )  =  <  XTOTAL  ( I )  *  <  K4-K6  )  + YTOTAL  ( I  )  * <  K1  *K6- K3*K4  )  + 

?,  ZTOTAL  <.!)*<  K3-K1.  )  )  /DET 

YBLUE  <  I  >  -  <  XTOTAL  ( I )  *  ( K5-K4  )  +YTOT AL  ( I )  >K  (  K2#K4-K  1  #K5  )  + 

S  ZTOTAL ( I ) *  <  K 1 -K2 ) ) /DET 

CONTINUE 

THIS  SECTION  FORMATS  AND  PRINTS  THE  COMPUTED  TIS I STIMULUS  0 A LUES  AND 

P  H  0  S  P  H  0  R  L  U  M I N  A  N  C  E  S 

IF  (STATUS  .EG*  ' DONE ' )  THEN 
WRITE (3, 810) 

WRITE <3, 820) 

DO  830  1=1 ,N 

WRITE  (  3  840 )  I  XTOTAL  <  I )  v  YTOTAL  <  I  )  ZTOTAL  ( I  ) 

CONTINUE 
WRITE < 3,050) 

WRITE (3, 360) 

WRITE (3, 870) 

DO  830  1=1 ,N 

W IS  I T  E  <  3  ,  8 4 0  )  I. Y IS E Ii  ( I )  Y G IS F E N  (  I )  ,•  YDI.  U! .  ( I  ) 

CONTINUE 
WRITE  <  3  >  850 ) 

WRITE (3, 850) 

END  IF 

FORMAT  (I  4X,  'THE  TIS I  STIMULUS  VALUES  ARE') 

FORMAT  (  2X ,  '  PO TNT  '  f  ?X ,  '  X '  ,  1.4X  y  '  Y '  ,  1  4X ,  '  7. '  ) 

FORMAT ( 3X , 12 , 7X , FB . 3  y 7X ,  18 . 3 , 7X  r F8 . 3 ) 

FORMAT ('  ') 


TOTAL ( I ) 
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Although  the  solution  is  now  optimal,  the  program  continues  until  the 
step-size  is  less  than  .0001  at  which  time  the  step-size  is  reset  back  to 
.5  and  point  number  three  is  fixed  in  its  place,  allowing  the  algorithm  to 
concentrate  on  trying  to  optimize  the  other  three  points.  The  step-size 
then  gradually  decreases  back  to  .0001  and  another  point  is  fixed,  and  so 
on,  until  all  four  points  are  fixed.  At  that  time,  the  locations  of  the 
points  are  stored  in  memory  and  all  the  points  are  unfixed  to  see  if  any 
more  improvements  can  be  made.  The  program  then  goes  through  the  whole 
process  again,  reducing  then  resetting  the  step-size  and  fixing  the  points 
until  all  four  points  have  been  fixed.  Because  no  improvements  could  be 
made,  the  program  is  done  at  that  time. 
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